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PREFACE 


Since  July  I,  1975,  Lawrence  Livermore 
I aboraloiy  (111)  has  been  participating  in  the  High 
Altitude  Pollution  Program  sponsored  b\  the  ll.S 
IX'partntent  ol  Lransportation’s  Federal  Aviation 
Administration  I his  report  describes  the  major 
accomplishments  and  significant  findings  during  the 
fiscal  year  ending  Septembei  JO,  1978,  for  work  per- 
tormed  at  III  under  Reimbursable  Agreement 
1 H)  I - L A7p\\  A 1-65 J I wo  major  research  areas  are 
covered  bv  this  agreement  ( I ) numerical  modeling  ot 
the  atmospheric  response  to  stratospheric  perturba- 
tions.  and  ( -I  the  processing,  archiving,  and  analysis 
ot  satellite  o/one  data  Progress  on  the  second 
research  area  has  recently  been  published  in  a report 
bv  I ovill  <•/  ill  ( I97g)  Only  work  performed  in  the 
lit st  area  will  be  reported  here  I he  successful  accom- 
plishment ol  the  main  subtasks  within  the  numerical 
modeling  area  has  required  contributions  and 
cooperation  from  many  participants.  1 he  work 
reported  here  is  the  result  of  the  collectiv  e effort  of  all 
those  listed  below. 

Scientific  Administration 

Joseph  B Knox.  Division  I eadet 
Frederick  M I uther.  Principal  Investigator 

Participants 

Julius  S.  ('hang 
William  H.  Ducwer 
Joyce  F.  Penner 
Raymond  I larp 
Donald  .1.  Wuebbles 
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ANNl’AL  RKPORT  OF  LAWRENCE  I 1YERMORE  LABORATORY  TO 
THE  FA  A ON  THE  HIGH  ALTITLOE  POLLl  TION  PROGRAM  - 1078 

1.  INTRODUCTION 


1 he  High  Altitude  Pollution  I’logiunt  til  MM') 
was  initiated  bv  I ho  f edeial  Aviation  Administration 
to  ensure  that  aircraft  engine  emissions  in  the  stratos- 
phere will  not  result  in  unacceptable  eflects  on  the 
biosphere  I awrcnce  I ivcrutoie  I aboratorv  till) 
has  participated  in  M M’l’  since  lulv  I'D.'  1 he  pri- 
mary research  emphasis  at  111  is  on  numerical 
modeling  ot  the  atmospheric  response  to  stiatos- 
phcnc  perturbations  A liiiulamental  tool  in  the  I I I 
effort  has  boon  the  one-dimensional  transport 
kinetics  model  We  are  also  m the  ptoces\  ot  de\  elop- 
ing a two-dimensional  transport  kinetics  model 
Suppoit  lot  model  development  is  shared 
between  HAIM’  and  anothei  protect  at  I 1 I involv- 
ing assessment  ot  the  chemical  and  climatic  eflects  ot 
atmospheric  nudeai  explosions  1 he  latter  protect  is 
tunded  bv  the  Division  ot  Mihtatv  Application 
(I'M  At  ot  the  Department  ot  l nergv  Although  the 
same  numerical  models  are  used  toi  both  studies,  the 
applications  ot  the  models  are  quite  ditterent  l or 
example,  the  study  ot  the  climatic  effect  ot  nucleai 
explosions  is  primaritv  concerned  with  the  time- 
dependent  response  ot  the  atmosphere  to  pulse  infec- 
tions ot  NO,,  whereas  the  II  A I*  I’  studv  is  concerned 
w ith  steadv  -state  inteciioi.s  ot  NO,  and  olliei  species 
In  August  I'J'b.  I l l "s  participation  in  II  AIM’ 
was  extended  to  include  a feasibilitv  studv  to  detei 
mine  whether  good  quality  total  o/one  values  could 
be  derived  from  infrared  measurements  bv  a multi- 
channel tiller  radiometer  (MHO  sensoi  carried 
aboard  a series  ot  satellites  operated  bv  the  I S An 
force  Infrared  radiance  data  became  available  m 
March  Id’?  trout  the  first  ot  a new  series  ot  satellites, 
which  are  a pait  ot  the  Detense  Meteorological 
Satellite  Program  Block  AD  series  I he  MI  R sensoi 
aboard  these  satellites  is  unique  in  that  it  is  the  lust 
cross-track  scanning  sensor  capable  ol  ptoviding 
o/one  measurements,  thus  increasing  the  amount  ot 
o/one  data  collected  I he  instrument  has  a lugliei 
resolution  than  previous  instruments. and  it  provides 
global  coverage  dav  and  night  Since  nud-ld",  this 
series  ol  satellites  has  provided  the  onlv  available 
global  o/one  data  I he  second  satellite  began  trails 
nutting  MI  R measurements  in  July  Id?',  and  the 
third  in  the  series  began  tiansmitting  MI  R measure- 
ments in  September  Id's  Ihe  fourth  (and  last) 
satellite  in  the  series  is  vet  to  be  launched  I he 


designed  lifetime  ol  each  satellite  is  two  years  I he 
o/one  data  derived  from  these  satellite  measure- 
ments aie  potentially  useful  to  the  scientific  com 
munity  tot  assessing  changes  in  the  o/one  lavei  that 
max  be  due  to  natural  oi  man-made  causes 

Ihe  feasibilitv  studv  was  completed  in  June 
Id's,  with  the  successful  ptocessmg  of  ?0  davs  ot 
data  taken  during  Id"  | he  qualitv  ol  the  data  was 
demonstrated  bv  comparison  with  corresponding 
o/one  data  obtained  at  selected  stations  m the  world 
sui  lace  netw  oi  k ol  Dobson  obsei  v atones  A descrip- 
tion ol  the  methodology  used  and  a presentation  ol 
the  results  ol  the  teasibihty  studv  aie  contained  in  a 
report  bv  lovill  «•/  ill  (Id'S)  Ihe  protect  is  now 
undctgotng  a transition  into  the  operational  phase 
during  which  all  ot  the  MI  R data  from  the  satellites 
w ill  be  piocessed  and  made  available  lor  distribution 
to  the  scientific  coniniumtv 

In  this  annual  report  we  describe  I I I \s  maiot 
accomplishments  and  findings  since  October  I.  Id" 

I arliei  H AH’  woik  is  described  in  om  previous 
annual  reports  (I  utlier  <•/  u/.,  Id'b.  |0”)  lu  the 
present  report  we  discuss  the  icsults  in  onlv  the 
atmospheric  modeling  area,  the  woik  on  satellite 
o/one-data  processing  and  analysis  has  recentlv  been 
reported  separately 

A modeling  effort  is  important  to  H AIM’  lot  two 
pimiarv  reasons  fust,  modeling  constitutes  a means 
ot  relating  engine  emissions  and  other  man-made 
perturbations  to  changes  in  atmospheric  composi- 
tion Because  ol  the  complex  physical  processes  and 
feedback  mechanisms  involved,  it  is  essentially 
impossible  to  assess  the  eflects  ol  emissions  without 
incorporating  these  processes  into  a model  and 
allowing  them  to  interact  Second,  modeling  contri- 
butes to  the  interpretation  and  analysis  of  experi- 
mental results  and  helps  determine  the  relative  prioi 
ities  lot  earning  out  this  woik 

Due  to  the  evei -increasing  complexity  of  the 
svstem  ol  governing  phy  sical  and  chemical  processes, 
theoretical  models  of  atmospheric  twice  species 
distributions  have  become  essential  tools  both  ding 
nosticallv  and  pt  ognosticallv  Such  models  are  based 
on  a set  ol  conservation  equations  describing  the 
cited  on  individual  trace  species  of  atmospheric 
tiansport,  chemical  and  photochemical  kinetics,  and 
various  natural  and  anthropogenic  sources. 


I 


* 


Since  the  detailed  modeling  ot  atmospheric 
chemist  rv  is  a relatively  voting  science  and  since 
overall  mlotmation  is  limited,  although  tapidlv 
increasing,  most  models  have  heen diagnostic  in  then 
scope  I he  most  comprehensive  validation  of  all 
classes  ol  models  lematns  diagnostic  W ithin  limits, 
the  distribution*  and  variabilities  ol  most  ot  the  tiace 
species  have  been  modeled  ivasonablv  well  * 

When  the  question  ol  antlu opogemc  intluence 
ot  the  o/one  lavct  was  i aisevl,  theie  was  a demand  loi 
the  models  in  prognostic  applications  1 he  piog 
nostic  capabilitv  ol  the  models  lot  making  piedic 
tions  over  decades  ot  centimes  has  tievet  been 
demonstiated  and  in  all  likelihood  will  nevei  be 
demonstrated  without  a well  developed  long  teim 
atmospheric  monitoring  progiam  ot  compaiable 
duiation  I muted  piognostic  validations,  such  as  the 
model  predictions  ol  the  ettccts  ot  I’olai  lap  \bsoi 
ption  t IV  \l  events  on  polai  sitatospheiic  o/one  and 
the  current  stratospheric  chloiotluoiomethanc  dis 
tnbutions.  are  reassuring  but  not  conclusive  In  pai 
ticular.  the  studv  ol  IV  \ events  tequttes  reevamina 
lion  because  ot  the  stiong  coupling  ot  Nik,  I'll',, 
and  HO,  cvcles  in  the  cuitent  chcmistiv.  which  was 
not  accounted  toi  in  the  onginal  studv  In  a sense 
t tv >t v studv  with  one-dimensional  models.  Ouewei  ei 
itl  tl')”l  demonstiated  that  undei  limited  ciicum 
stances  diagnosticallv  similat  one  dimensional  mod 
els  ol  the  stiatospherv  can  indeed  vield  veiv  dilleient 
predictions  on  the  etlect  ot  NO,  ettluents  tiom  SS  I N 
on  o/one  I'onsequentlv . it  is  tuvessatv  to  tecogni/e 
the  distinction  between  piognostic  and  diagnostic 
applications  ol  the  same  model  and  allitm  then 
individual  limitations  and  uncertainties 

Oiagnosticallv  and  prognosticate.  box  models 
can  ptovule  only  limited  mloimation  In  fact,  vine  to 
its  simple  structure,  m certain  applications  it  is  lughlv 
sensitive  to  uncertainties  m input  variables  I 01 
example,  analvsts  ol  atmospheric  lifetime  ol  long 
lived  halocarbons  with  box  models  is  lughlv  sensitive 
to  uncertainties  m the  ratio  ot  measured  global  ahun 
dance  to  the  total  amount  released  I ot  I 1 1 and  I 
1 2,  an  uncer  taint  v ot  V ItV  , m this  latio  w on  Id  lead  to 
an  uncertain! > ot  25-50*7  m the  deducevl  atmosphei  to 
lifetime  (Chang  and  I’enttet.  I*'"$!  I'ndeistandmg 
this  type  ol  uncertainty  is  ot  particular  importance 
tor  the  box  model  because  ot  its  application  m 
lifetime  analvsts  Similar  uncertainties  aie  in  the 
more  complex  models  but  then  applications  are 
normally  hevond  such  simple  analvsts  \t  piesent 


• 1 he  mo5|  noUiMe  exception*  nuf  K the  hi||h  l‘IO  meaMitcmcnl* 
\n*1n  *on  \ 


theie  exist  seveial  usetul  techniques  toi  the  analvsts 
ol  box  model  sensitivity  and  uncertainty  Ihcteloie 
one  would  expect  such  an  analysis  pilot  to  the  mtei 
piotatiou  <>l  modeling  icsults  Hox  models,  with  oi 
w ithout  associated  sensitiv  itv  ana Iv  sis.  can  best  set  \ e 
to  establish  the  telative  importance  ot  phvsica! 
processes  and  suggest  those  pioblems  that  ate  most 
woitliv  ot  detailed  examination  with  moie  complex 
models 

I he  stiuctute  ol  one -dimensional  models  has 
Iveen  analv/ed  in  detail  m the  lepoit  In  the  National 
Research  I'ouncil  y l '’f>b>  and  in  the  NASA  ivpoit 
(Hudson,  |d”t  I he  veitical  tianspoit  parameten 
ration  is  based  on  globallv  aveiaged  values  ot  veitical 
disii  ihutions  ot  tiacets  such  as  N O and  i'll. 
I'onsequentlv.  lot  othei  piedommantlv  upwaul 
moving  tiacets  the  one  dimensional  model  should 
Iv  considered  to  be  a globallv  aveiaged  model 
this  is  piobablv  also  tine  loi  the  net  downward 
moving  ttaceis  m the  altitude  i.ingc  ol  the  tiopo 
pause  Hut  m the  siiatospheie  the  detailed  icpie 
sentation  ot  the  photochemical  mteiactions  is 
such  that  onlv  local  vauables  such  as  tempciatuie 
and  solai  llux  mtensitv  (both  in  space  and  timel  can 
be  used  1 his  is  a necessaiv  condition  that  sigmti 
cantlv  complicates  the  mtei  pi  etatn'ir  ot  model 
icsults  lire  local  chemical  litetime  ot  all  stiatos 
phene  chemical  tiace  species  depends  on  the  Uval 
cv'nccntrations  ot  o/one  and  the  closelv  ev'iipled 
atomic  oxvgens  iH  I't  and  l'(  1"  1 Irese  concentia 
tions  aie  duectlv  piopomonal  to  the  local  solai  this 
As  a result,  the  sell  cousistentlv  derived  oikl  oxvgen 
concentrations  m the  one  dimensional  model  aie 
positivelv  correlated  with  the  seasonal  solai  this 
mtensitv.  eontiatv  to  obseivation  ot  the  seasonal 
vauation  ot  hemisphencaland  global  o/one  distnbu 
tions  In  the  absence  ot  detailed  mloimation  on  the 
spatial  and  temporal  vauations  ot  all  the  tiace 
species,  eithei  tiom  measuiements  oi  mote  complex 
models,  it  is  not  possible  to  detive  the  piopei  avei 
aging  (weighting!  taclois  toi  the  nonlmeai  photo 
chemical  and  chemical  interaction  ptocesses  in  the 
one  dimensional  model  1 heiefote.  all  one  dimen 
sional  models  ate  limited  to  ivpiesentmg  local  siia 
tospheric  conditions,  nr  particular  that  latitudinal 
region  and  season  where  the  model  domed  o/one 
profile  closelv  repivsents  the  Uval  measuiements 
lvpicallv.  one  dimensional  models  host  ivpiesent 
the  photochemist  tv  ot  about  hl'N  latitude  at  about 
the  tall  equinox  1 hat  this  is  not  completelv  satistac- 
toiv  is  quite  obvious  since  this  is  also  the  region 
w here  there  i'  consideiahle  stiatospheiic and  tiopos 
phene  exchange,  hens'e  we  exjvct  local  coiurntia 


transport  processes.  Consequently,  because  of  this 
mix  of  partially  local  and  global  representations, 
both  the  diagnostic  and  prognostic  application  of  the 
one-dimensional  model  must  be  applied  and 
analyzed  with  great  care  with  regard  to  the  intrinsic 
compatibility  of  the  model  and  its  intended  use,  case 
by  case. 

The  results  of  a model  calculation  can  be  evalu- 
ated by  comparison  with  suitably  obtained  atmos- 
pheric measurements.  Several  basic  types  of 
information  have  proven  to  be  useful:  height  and 
latitudinal  distributions  of  individual  species  and 
groups  of  species,  partitioning  of  related  species,  and 
local  diurnal  and  seasonal  variations  of  individual 
species.  These  comparisons  when  available  consti- 
tute a necessary  test  mostly  for  the  diagnostic  end  at 
times  prognostic  adequacy  of  the  model.  Interpreta- 
tion of  the  results  from  such  tests  are  not  always 


simple.  Many  of  the  comparisons  must  be  accepted 
with  degrees  of  qualification  It  is  the  overall  rea- 
sonableness of  the  models  and  their  predictions 
that  provide  a high  level  of  confidence  in  their  basic 
correctness.  Of  course,  this  is  often  a matter  of 
judgment. 

T he  measured  atmospheric  trace  species  con- 
centrations. which  constitute  the  data  base  for 
comparison  with  model  results,  have  shown  a 
significant  level  of  local  variability.  This  atmospheric 
variability  has  caused  such  comparisons  between 
measurement  data  and  modeling  results  to  be 
influenced  by  judgment  and  at  times  provided  a 
ready  excuse  for  inconclusiveness.  As  data  base  and 
theoretical  models  continue  to  improve,  a more 
precise  formulation  of  atmospheric  variability  and 
its  influence  on  data  interpretation  must  be 
developed  and  applied. 


2.  MODEL  SIMULATION  OF  OBSERVED 
STRATOSPHERIC  VARIATIONS 


2.1  The  Natural  Stratosphere 

A fundamental  tool  in  the  LLL  effort  has 
been  the  one-dimensional  transport-kinetics  model 
(Chang  and  Wuebbles,  1977;  Luther  el  al..  1977). 
This  model,  which  includes  as  complete  a set  of  the 
important  chemical  and  photochemical  reactions  as 
is  feasible,  is  designed  for  time-dependent  perturba- 
tion and  sensitivity  studies.  The  model  i icludes  30 
chemical  species  and  97  chemical  and  photochem- 
ical reactions.  The  chemical  rate  coefficients  cur- 
rently used  in  the  model  are  listed  in  Table  A-l 
(Appendix  A).  Stratospheric  water  vapor  concen- 
trations are  now  computed  at  altitudes  above  12  km 
in  our  standard  model,  although  a fixed  concentra- 
tion profile  may  be  specified  as  a user  option. 

Species  concentrations  are  computed  at  44 
levels  in  the  atmosphere,  extending  from  the  ground 
to  an  altitude  of  55  km.  The  model  uses  an  accurate 
numerical  method  for  solving  stiff  systems  of  differ- 
ential equations.  Vertical  transport  is  parameter- 
ized using  a one-dimensional  diffusion  formulation 
that  describes  hemispheric-average  net  vertical 
transport  by  an  altitude-dependent  diffusion  coeffi- 
cient. The  Chang  (1976)  diffusion  coefficient  profile 
was  used  for  all  of  the  calculations  presented  in  this 
annual  report.  The  model  can  include  temperature 
coupling  between  changes  in  composition  and  reac- 
tion rate  coefficients,  but  this  feature  is  not  used  in 


the  standard  model.  The  standard  model  includes  a 
totally  self-consistent  diurnal  averaging  procedure 
for  nonlinear  chemical  reaction  terms  (Luther  et  al.. 
1977).  Unless  otherwise  stated,  the  standard  model 
was  used  for  the  calculations  presented  in  the  follow- 
ing sections.  Changes  made  to  the  model  during  the 
past  year  are  described  in  Appendix  A.  In  this  section 
we  describe  the  natural  stratosphere  as  computed  by 
the  model.  The  sensitivity  of  the  model  to  various 
stratospheric  perturbations  is  described  in  Sec.  3.1. 

The  ambient  stratospheric  species  concentra- 
tions that  are  presented  here  correspond  to  a simula- 
tion of  present-day  conditions.  The  historical  release 
of  CFM's  is  included  assuming  constant  production 
since  1975.  Since  the  chlorine  budget  of  the  stratos- 
phere is  continuously  increasing,  this  time-depen- 
dent calculation  not  only  provides  a current  stratos- 
phere for  diagnostic  purposes,  but  it  also  provides  a 
test  for  the  model  in  predicting  the  accumulation  of 
CFM’s  in  the  troposphere  and  their  vertical  distribu- 
tion in  the  stratosphere.  In  this  calculation  it  is 
assumed  that  presently  observed  tropospheric  con- 
centrations of  CCIj  and  CHiCI  are  the  sources  for 
background  chlorine.  The  species  CH>CCL  is  not 
used  because  of  difficulties  with  data  precision  and 
accuracy  and  the  apparent  conflict  with  measured 
HO  concentrations  in  the  lower  troposphere  (see  the 
discussion  on  CHiCCli).  With  or  without  CHiCCL 
as  a source  of  present  stratospheric  chlorine,  the 
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analysis  ol  the  tollowmg  sections  is  the  same  with 
certain  exceptions  a>  noted. 

1 he  computed  stratospheric  distributions  ol 
minor  trace  species  corresponding  to  .human  I *>78 
are  shown  in  Figs  I through  b.  I hese arealldiurnally 
(24-hr)  aseraged  quantities  l he  relevant  noontime 
and  nighttime  values  are  discussed  in  individual  sec- 
tions I he  performance  ol  stratospheric  models  in 
picdicting  the  vertical  distribution  ol  known  trace 
species  that  are  believed  to  be  important  to  o/one 
chemistry  are  discussed  in  the  following  sections 
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M(i.  I Sprite  cont'rnlralitin  |Hi>rilr>  for  O «tul  <><  IMiomiuilrti 
f«*  Jtnuan  I*?*. 


I race  Species  With  l emu  Stratospheric  C hemical 
lifetimes 

The  trace  species  N.O,  C'F.C'I;,  CFC'h,  t'C'h, 
t Hit  ('h.  ('Hi.  I'HiCI,  and  H.-O  all  have  stratos- 
pheric chemical  lifetimes  ol  at  least  a year  Conse- 
quently. then  stratospheric  abundance  is  totally 
determined  bv  (luxes  Irom  the  troposphere  (strictly 
speaking  this  is  not  true  for  H O as  will  be 
explained)  In  tact,  these  species  are  the  sources  for 
all  the  mmoi  trace  species  in  the  stratosphere  that 
are  known  to  be  essential  in  o/one  chemistry 

N20 

I he  basic  input  variable  for  modeling  N.-O  is  its 
observed  tropospheric  mixing  ratio.  Models  can 
either  use  this  number  directly  or  convert  it  to  a 
source  llux  derived  Irom  the  measured  mixing 
ratio  and  an  assumed  atmospheric  lifetime.  Asa 
result,  agreement  with  observations  in  the 
troposphere  is  expected  (.'20  ppbv  vs  .'25  * 
10  ppbv).  In  the  stratosphere,  the  N;()  local 
chemical  lifetime  increases  from  b months  at 
.'5  km  to  thousands  ol  years  near  the  tropo- 
pause  In  fact,  one  may  state  that  the  variability 
ol  its  vertical  concentration  profile  in  the  stra- 
tospheie  is  strictly  due  to  atmospheric  transport 
processes.  One-dimensional,  and  even  two- 
dimensional  models,  can  at  best  prov  ide  a global 


M(<  ' I race  cpectes  mixing  ratios  computed  (or  Januart  I*’* 


qualitative  picture  hut  no  exact  quantitative 
agreement  locally  As  is  well  known,  one- 
dimensional  models  can  indeed  simulate  the 
average  veitical  distribution  ol  N O It  should 
be  emphasized  that  both  N O and  t'H4  have 
been  used  as  tracers  tor  deriving  the  etlective  net 
vertical  dittusion  coelticient  in  one-dimensional 
models.  I hcreJore,  such  agreement  between 
data  and  model  results  is  a necessary  condition 
It  confirms  that  the  model  is  functioning  as 
expected,  i e . the  more  complex  lull  model  does 
not  change  the  expected  chemical  cycles  of  VO 
and  t il.,  hence  no  reevaluation  of  the  effective 
diffusion  coefficient  is  necessary 


A.  theoretical  and  measured  concentrations  of  Oil 
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CH4: 

Methane  behaves  much  the  same  as  N.-O  in  the 
atmosphere,  except  that  it  has  a comparatively 
short  lifetime  in  the  troposphere.  The  model 
used  1.31  ppmv  based  on  earlier  data,  which  is 
20%  less  than  current  data  (1.65  ± 0.1  ppmv). 
This  leads  to  similar  underestimates  throughout 
the  stratosphere.  It  also  leads  to  lesser  HO,  from 
(Tl4,  which  would  have  a small  effect  on  the 
effectiveness  of  added  chlorine  on  ozone.  Since 
the  local  chemical  lifetime  of  CH4  increases 
from  6 months  at  40  km  to  more  than  10  years  at 
the  tropopause,  its  stratospheric  variability  is 
also  controlled  by  transport  processes.  In  par- 
ticular the  qualitative  latitudinal  similarity  in 
the  vertical  profiles  of  N.O  and  CH4  from  5°S 
and  40°  N confirms  this  view. 

CCI4  and  CH3CI: 

The  model  used  approximately  the  measured 
concentrations  as  input  parameters.  The  tropos- 
pheric concentration  of  CCI4  is  1 12  ppbv.  which 
isat  least  10%  lower  than  the  currently  available 
averages  but  certainly  within  the  limits  of  data 
uncertainty.  The  CCI4  is  not  destroyed  in  the 
troposphere  but  begins  to  photodissociate 
quickly  in  the  stratosphere.  Its  local  chemical 
lifetime  is  of  the  order  of  6 months  at  about 
24  km.  Therefore,  in  addition  to  the  instrumen- 
tation difficulties  associated  with  grab  samp- 
ling, there  should  be  a high  level  of  transport 
variability  which  makes  stratospheric  measure- 
ment of  CC’I4  and  its  data  interpretation  a rather 
difficult  task.  The  CHiCI  measurement  data 
have  shown  such  fluctuations  that  it  is  difficult 
to  assess  the  adequacy  of  model  results.  Never- 
theless. with  an  input  of  746  pptv  at  the  surface 
(which  is  approximately  in  the  middle  of  the 
range  for  Rasmussen's  data  in  1977),  the  model 
estimates  580  pptv  at  the  tropopause.  which 
seems  to  be  consistent  with  the  550  pptv  re- 
ported by  t'ronn  <7  al.  (1977).  The  estimated 
local  chemical  lifetime  is  approximately  I year 
up  to  26  km.  This  again  makes  transport  varia- 
bility an  important  factor  in  the  interpretation 
of  ('ll  i('l  data. 

CF2CI2  and  CFCI3: 

The  total  abundance  and  vertical  distributions 
of  F-l  I and  F-12  serve  as  true  tests  of  the  ability 
of  simple  one-dimensional  models  to  represent 
the  average  net  vertical  transport  of  trace 
species.  Together  with  CHiCCii  the  release 
rates  of  these  man-made  chemicals  arc  well 
documented.  Since  the  model  uses  the  time- 
dependent  release  rate  in  computing  its  atmos- 
pheric distribution,  a comparison  of  both  the 


model-derived  tropospheric  and  stratospheric 
budgets  with  measurement  data  would  be  a 
stringent  test  of  the  adequacy  of  the  model. 
Global  averages  of  F-l  I and  F-12  in  the  tropos- 
phere are  not  easily  obtainable.  If  the  observed 
northern-hemispheric  averages  in  January  1977 
are  used  to  estimate  the  global  average  in  Janu- 
ary 1978,  then  the  model  is  at  most  10%  lower 
in  the  tropospheric  burdens  of  F-l  I and  F-12. 
About  5-8%  of  this  is  attributable  to  the  accu- 
racy of  the  source  function  approximation  used 
in  the  model.  This  then  provides  a partial  (but 
important)  confirmation  that  the  transport 
approximations  based  primarily  on  NT)  and 
CH4  indeed  have  a more  general  validity  and 
applicability.  Similar  to  other  trace  species  with 
long  local  chemical  lifetimes,  both  F-l  I and 
F-12  are  sensitive  to  variations  in  local  transport 
processes  (at  altitudes  up  to  30  km  for  CF'C’li 
and  up  to  40  km  for  CFjCIj).  This  is  again  con- 
sistent with  observation.  In  general  the  model- 
predicted  stratospheric  vertical  profiles  of 
CF'.'C'b  and  CFCIi  are  both  slightly  on  the  high 
side  of  the  data  averages  (taking  the  model  tro- 
popause to  be  at  12  km).  The  model  seems  to 
show  slightly  more  destruction  in  the  lower  stra- 
tosphere (5-10  km  above  tropopause)  forCFCIi 
than  the  data  suggest.  On  the  other  hand,  this 
could  be  due  to  weaker  net  upward  transport  in 
the  model.  Fortunately,  in  the  regions  of  maxi- 
mum destruction  of  CFC'h  and  CF.-Cb  (20-30 
km  and  25-35  km,  respectively),  the  model- 
predicted  profiles  are  in  reasonable  agreement 
with  available  data. 


CH3CCI3: 

At  this  time  there  is  considerable  uncertainty 
in  the  interpretation  of  data  and  the  modeling 
of  atmospheric  concentrations  of  CHit'Cli.  A 
simple  box  model  analysis  based  on  measure- 
ment data  would  yield  a CH1CCI1  atmospheric 
lifetime  of  5-12  years,  which  in  turn  would  imply 
an  average  tropospheric  HO  concentration  of 
less  than  10"  molecules  cm1.  However,  most 


coupled  stratospheric-tropospheric  models 
have  more  HO  than  this  estimate,  which  then 
leads  to  a much  lower  C'HiC'Cl,  concentration 
in  the  atmosphere  if  the  well-known  release 
rates  are  used  (Neely  and  I’lonka,  1978).  Fur- 
thermore. there  is  no  generally  accepted  set  of 
HO  measurements  in  the  troposphere.  Al- 
though CHiCCii  is  estimated  to  be  a major 
source  of  stratospheric  chlorine  and  could  pro- 
vide another  direct  test  on  the  adequacy  of 
models  at  this  time,  it  is  not  constructive  to 
include  it  in  the  one-dimensional  model  in  the 
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same  inannci  as  the  othei  source*  ot  stratus 
phene  chlorine  I lie  model  vise'll  foi  these  ealeti 
lations  el ui  not  mcluele  this  souree.  which  eonUI 
eontuhiite  from  0.2  toO  I ppbs  ol  ehlonne  line 
omission  is  ol  lie'  practical  consequence  te>  the 
present  Minis 

11)0 

Methane  is  the  pnncipal  siiurceol  siraiosphenc 
watei  At  the  iropopausc  II  t)  is  maintaineel 
both  b*  a dossnward  this  trotn  the  stratosphere 
ami  an  upwarel  Hus  Iroin  the  troposphere  At 
this  altitvule  the  tropospheric  stratospheric  ex 
change  process  seises  to  maintain  II  •()  at  the 
satutatiesn  level  I Ins  lixeei  beunnlars  coiulition 
helps  maintain  the  II  li  coiu'cntration  in  the 
iest  ol  the  stratospheie  With  such  control  it  is 
not  surprising  to  see  one  dimensional  models 
geneiatc  reasonable  li  t)  vertical  piotilcs,  siait 
mg  ssith  t ppms  at  the  tiopopause  and  slow  Is 
increasing  to  * ppms  at  '0  km  I his  is  ssell 
esithnt  the  iincerlainiv  ol  the  data  Again,  since 
stratospheue  ssatet  sapoi  has  onls  weak  chem 
teal  souiees  and  sinks,  its  sanahilils  is  mostls 
due  to  transport  lluetuations 
• 

I able  I suminaii/es  the  so  vines  ol  ehlonne  in 
the  moilel  and  compaics  them  ssith  estimates  tiom 
lioposphene  data  Vs  ssas  abends  pointed  out.  the 
model  underestimates  the  total  budget  bs  not 
including  the  uncertain  contribution  ol  methsl 
chlotoloim  V\  it h a total  I It),  mixing  latio  ol  ' ppbs 
in  the  tioposphere.  the  l 'It),  mixing  ratio  is  leduced 
to  I ' ppbs  at  40  km  due  to  the  natural  time  delas  in 
tipssatd  nanspoit  and  the  photochemical  destine 
non  ol  the  source  molecules  at  lowci  altitudes  In 
evunpaiison.  the  total  Nti,  mixing  ratio  at  40  km  at 
I s o ppbs . and  the  total  lit',  mixing  latio  at  40  km  is 
' ’ ppbs 


I able  I Source*  of  stratospheric  chlorine 
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l'race  Species  with  Intermediate  Stratospheric 
Chemical  I ifetimes 

I race  species  such  as  UNO,  and  IK'I  have  stia 
tospheue  chemical  lifetimes  ot  about  a week 
Similarly,  the  trace  species  families  O,  (O('D)  t 
Ot'lM  » Odaiul  NO, (N  t NO  t NO  t NO,  t UNO, 
t ’N  t),)  also  base  net  stratospheric  chemical  hie 
nines  ol  about  a week  I he  ehlonne  Iannis  CIO,  (Cl 
♦ CIO  t IICI)  has  a somew  hat  shoitei  stratospheric 
chemical  lifetime  on  the  order  ol  I ^ to  I das  I his  is 
due  to  the  strength  ol  the  coupling  reaction  CIO  * 
NO-  ♦ M • CIONtV  t M.  since  none  ol  the  abose 
(usual)  t atm  lie*  includes  CltlNt)-  II  sseconsulei  the 
broadei  Iannis  CIO,  (CIO,  t CIONtM.  (hen  it 
becomes  a pine  traeei  ssith  no  chemical  loss,  except 
possibls  thiough  the  set  to-be  completely  quantified 
icactions  involving  IIOt'1  In  olhci  words,  it 
simultaneous  measuiements  aie  made  ol  all  (he 
membeis  ol  the!  It ),  lannlv . then  all  vanahiliis  in  the 
data  must  be  vine  to  pits  steal  processes  othei  than 
photochemisiis  On  the  othei  hand,  il  onls  CIO, 
Inmtls  membeis  ate  measured  sunultnneousls  lie. 
without  CltSNO  ).  (lien  the  variations  due  to  chem 
icnl  leactions  involving  l IONO  must  beconsivleied 
I ogethei  with  IIO,  ( li  t HO  t II  t 'ICO-),  the 
t ).,  \t  s,  and  l It',  families  include  almost  all  (he  dace 
species  that  aie  ol  maiot  interest  in  stiatosphciic 
chemist i\  then  shod  net  sdatosphene  cheiiueal 
hleiimcs  signitv  the  dominance  ot  chenustiv  in  then 
hie  cvcles  and  thru  pnmais  siiatosphencoi  igm  I he 
(),.  list',  and  IICI  downs* aid  Duxes  (rom  the  siia 
tv'sphei e i epresent  almost  all  the  nanspoit  loss  ol  live 
dace  species  m then  respective  families,  and  the 
nanspoit  lime  scales  are  all  ol  the  ocdci  ot  a seat  oi 
moie  Because  ot  these  similarities  and  leladotts,  it  is 
appn'pnate  to  vliseuss  these  tatmhes  Iveie  and  leave 
live  paid!  lotting  vv  it  Inn  mdiv  idual  (ami  lies  io  (he  ncvi 
sections 

Os 

I oi  the  most  pail,  the  sdatosphene  O,  species 
ate  m the  lot  in  of  O,  I he  ( >(  'l’l  begins  to  con 
tubule  moie  than  a lew  percent  onlv  at  altitudes 
above  appioximatelv  4'  km  V comparison  ol 
model  domed  t distnbution  vv  it  Is  ntmosphei  ic 
■lata  actualls  icpreseuts  an  evaluadon  ol  the 
whole  O,  tanulv  m (he  snatospheic  I he  model 
denved  vcidcal  piolile  ol  t),  (I  ig  II  compares 
well  with  live  aseiage  profile  obtained  In 
VV  atanabe  and  I olunaisni  1 4 'Mat  hi  N llclovv 
's  km  the  O,  coitcouiiadoiis  ate  levs  than  the 
midlatitude  piolile  ol  Kiuegci  and  Mm/net 
(|4*(,|  | Ins  is  an  expected  icsvill  since  the  mid 
latitude  piolile  included  eondibutions  liom 
high  latitude  o/one  distributions  that  aie  the 
lesuli  ol  meridional  nanspoit  VV  ithm  the  one 
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dimensional  model  tiamewoik.  this  ttauxlatcs 
to  a low  1*1  aveiage  veilieal  dittiiMon  coclltcieiu 
m i he  tange  ol  I1*  20  km 

I he  piincipal  souice  ol  odd  mitogen  unite  Mia 
tospheie  is  VO  thiough  the  teaetnm  with 
()(  l>>  I here  in  also  a small  amount  piodmed 
Irom  cosmic  lav  honihaidnieni  St  nee  most  ot 
the  ptodnction  in  heh'w  in  km  with  a Ntttall  net 
Ii'nn  1 1 vim  the  \ ♦ NO  icaction  above  4'  km. 
total  NO.  in  taulv  mill  ninth  vliNtiihtilevI  between 
U>  4n  kttt  with  a mixing  tatio  v'l  1 4 pphat  40  km 
I he  cottcxponding  Ml'.  mixing  ratio  ol  10  4 
pph\  at  2'  km  eompaieN  well  with  ohNeivaiioitN. 
wheieai  at  " km  the  NO.  mixing  latio  liv'iti 
the  ntovlel  (I.V4  pplnl  in  ai  the  high  eml  v>t  the 
NiiatONpheue  data  lattge 

I hiv'iighoul  the  NitaloNpheie.  the  loeal 
elteimeal  hletimcN  ol  NO.  aie  v'l  the  oulei  v'l 
social  daw  l'itn.xcv|ucutl\,  all  (except  a nv-i  v'l 
exaeth  Nimullaueoiis'  meaxuientenlN  vit  NO. 
NO  ami  UNO*  eoiilvl  he  stihieel  lo  unknown 
mtlueneex  vine  to  elteimeal  aitvl  tiaiiNpv'ii  pivi 
cesses  I he  UNO,  comeiilialions  wunptilevl  h\ 
the  model  .ne  on  the  high  enle  v'l  iltevl.ua  i.itige. 
hut  tine  in  vpnle  tvpical  v'l  all  one  vhmeitNionnl 
mv'vleh  I he  UNO,  mixing  lain'  peakN  al  about 
' 28  km  ami  in  eoitNiNtent  w ii It  all  the  available 
vlata  I he  II  NO,  Ii.in  a loeal  eltemieal  hletnne  v'l 
10-20  daw  thutiiglivHii  the  altitudes  whoie 
II N 1 > . in  ntv'Nt  ahuitvl.ini  aitvl  ean  x ;n x sigmli 
eanth  v>\ei  a lew  vlaxx,  asvvas  Iv'tiitvl  h\  Muieiav 
<7  .//  ( |4*S)  Ouahtativeh.  the  tneteaxe  li'w.iivl 
the  pv'le  in  ei'ltimn  UNO.  in  eoiiNiNtent  with  the 
present  phv'toeltemieal  tlteorv.  Inn  iwv'  vlunen 
Niv'ttal  ntv'vlelN  have  itv't  pioxivlevt  Ntillteiemh 
ilelatlevl  v|iiantilative  .ih.iIvnin  \ei 

1 

Stitee  the  tv'ial  Nltatv'Npheite  chlorine  in  the 
ntv'vlel  ean  onlv  be  removed  aN  lli'l  ihiv'iiglt 
vlv'wuwaivl  tiaiiNpv'ii  aeiv'NN  the  ttopopaiixe.  un 
Ntiatv'Nphene  mixing  lain'  in  appi.'ximaielv  the 
vame  aN  that  v'l  the  souiee  NpeeieN  II  we  leieel 
the  three  high  valnex  ot  i'll'  veilieal  piotileN  h\ 

\ itvlvi nv'ii  tIO’Sh  then  the  tv'ial  i'll',  mixing 
ratio  ai  .'0  km  in  appioxintateh  <>  > 10  pphv. 
aitvl  al  O km  it  in  I h 2 ' pphv  I he  I I I ntv'vlel 
gives  0 'x  pphv  at  20  km  aitvl  I h'  pphv  at  " km. 
vvltieh  .tie  vpiite  ev'inpaiahle  to  the  expenmental 
values  I he  high  l it'  ptolilex  hv  thetusebv'x 
wv'iilil  imheate  a l It',  mixing  tali"  v'l  I S o 0 
pphv  v'l  ehlv'ime  al  ahoul  40  km  I ogelltci  with 
the  eNtuuaievI  10  1'  pphv  v'l  Mi  l.  this  wv'tilvl 
iev|iiue  a sv'tnee  ol  ' 10  pphv  v'l  tiv'poxphetit 
ehlv'ime  I Iiin  in  Nimplv  itv't  aeev'iintevl  lv«i  hv 


eiiiient  umlei Niamling  ol  the  wtuieev  aitvl  viukv 
v'l  ehlv'itne  in  the  atmoNpheie  \lt  hough  \V  aieiv 
.7  ,it  1 1 4 ’S'  pii'v  ivlv  il  Nvutie  vlata  xiiggeNiing  vIin 
agieentent  with  vuie  >'l  the  ilnee  high  piv'lile 
uieaNiuemeniN.  a ntoie  vlelimlive  explanaitv'it 
v'lt  the  Nv'tuee  v'l  eiiv'i  eilltei  v'lt  the  photoehem 
leal  lliei'iv  v'i  the  itieaNUieinentN  hv  \nvlei noii  in 
Mill  iev|uneil  In  the  h'wei  Miatospbeie.  lit  I o 
the  pievlv'itnnani  Iv'i in  v'l  l It '.  I hv  ntv'vlel  givev 
t'  2S  pphv  at  I'  kin  aitvl  0 " pphv  at  20  km 
vvltteh  ean  he  w'lttpaievl  i>'  the  meaNtuevI  valuev 
v't  0 2 pphv  aitvl  0 I 0 ’ pphv  at  I'  aitvl  20  kin. 
ivepeetivelv  II  tlieve  valuev  aie  iivevl  h'  estimate 
the  vlv'wuwaivl  tin \ v'l  v'vlvl  ehli'iuie.  then  thev'iv 
aitvl  expel  intent  aie  within  a laetot  v'l  ' ' I Iiin 
again  piv'v  nlev  anoiltei  utvlepeitvlent  estimate  v'lt 
the  uneenamiv  vine  iv>  liaiiNpoit  I lie  dip  v'l  the 
111  I veilieal  piv'lile  al  '0  km  as  tv-potted  hv 
I a/iue  v’i  >il  1 1"  "I  aitvl  \\  tlliatUN  , i ,tl  1 1"  '(>'  i' 
nv'i  explainable  vviihut  the  ev'ttlexi  v'l  one 
vliiueuNii'tial  niv'vlelv  Mui  a pifhiuuiaiv  iwo 
vluneitNn'ital  > alv  nlativ'ti  hv  I’v  le  1 14  'si  suggests 
that  III  1 has  loeal  maxima  aitvl  minima  m (lie 
v eilieal  d u eel i v'lt  extettvlutg  liv'iti  ueai  the  evpia 
Iv'i  lo  the  polai  legion  II  v'lte  vveie  i »'  llv  a 
hallv'v'it  aiv'iuivl  <2  N.  u would  ipme  likelv 
sample  iluough  this  minimum  legn'it  aitvl 
obtain  a dip  in  the  veilieal  piv'lile  I Iiin  ieMilt 
ev'iiplv-vl  with  the  t'lO  tesull  (see  next  section' 
luithet  vlentv'itNiiaies  the  limitattiuts  v'l  a v>ne 
duneitstv'ital  ntovlel  aitvl  ai  the  same  tune  ex 
pi  esses  a stioug  vole  v'l  confidence  Iv'i  the 
potential  heuelu  ol  iwv>  vlnueitNional  models 


1 1 mcc  Species  with  Sltoit  Stmtospherie  l hemivnl 
I del  tines 

I he  ItavX'  Npv'en-N  \t',  NO-,  Nl'v,  N v'.,  i'|  I'l. 
I 'l.  I'll'.  I'li'NO  III'.  Hi'  aitvl  Hi'  all  have  much 
sltv'itei  eltemieal  lilettmes  than  the  piesioiislv  »l in 
cussed  species  \s  a lesuh  i Iteu  utvliv  ivlual  Iv'i  a I con 
eentiaiiv'its  .tie  vpme  vauahle.  nstiallv  v'n  the  tune 
scale  v'l  It. 'ins  lo  seconds  I hev  all  vaiv  dttituallv 
(except  h'i  II  t'  season. tils  and  piobablv  eveitv'n 
the  tune  scale  v'l  ehaugiug  alhcvlv's  vim’  iv<  the 
piescme  v'l  aeii'sv'ls  ami  i h'iiits  at  h’w  ei  almiivles  Vs 
is  ex|'v  % led.  ii  is  not  possible  (ami  nstiallv  impossible' 
to  alvvavs  obtain  local  conditions  such  .t'  solai  this 
Iv'ittpei  aim  v aitvl  local  com  cult  at  ions  v'l  othei  u a.  e 
'peciv'N  when  doing  a partieulai  measuicmenl 
II.' ii. v.  it  is  taiv'lv  possible  I.’ ijuaniilauv civ  inlei piei 
the  av  a liable  data  on  these  nave  spe.  ies  vv  ult  iltevMv 
tieal  ntv'vlel'  Hv'wevei.  il  individual  collections  v'l 
vlata  aie  used  as  a gioup  then  I Its-  total  local  v ai  i.tbil 
uv  in  the  vlaia  in  v'Heu  a gv'v'd  guivling  viilein'it  Iv'i 
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testing  the  theory  I uithetmore,  variability  m t he 
shape  nl  individual  \ ertical  profiles  is  also  a uselul 
iiulieatoi  ol  the  domain  ol  pin  steal  results 
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NO.  NOj,  NO,  and  N2Os 

Since  1972  there  have  been  a considerable 
number  of  NO  measurements  in  the  strato- 
sphere. By  now  its  stratospheric  variability  is 
well  established  and  any  attempt  to  model  a 
particular  set  ol  data  has  long  been  deemed  as 
relatively  unimportant,  except  tor  the  diurnal 
variations  of  the  NO.  species  As  was  pointed 
out  m our  previous  work  (1  uther  <7  <;/ , l*)7<>), 
successful  modeling  of  the  NO,  diurnal  varia- 
tions confirmed  locallv  some  of  the  chemical 
mechanisms  important  in  the  catalytic  cycles 
affecting  o/onc  Our  one  dimensional  model 
yields  an  NO  mixing  ratio  of  0 I ppbv  at  20 
km.  graduallv  increasing  to  *)  0 ppbv  at  >0  km 
I Ins  theoretical  profile  is  approximately  m the 
middle  of  the  available  data  Without  a 
simultaneous  measurement  of  NO-,  individual 
NO  data  will  be  of  steadily  declining  use 
fulness  m the  future 

Modeling  ol  stratospheric  NO-  distnhu- 
tions  suffers  the  same  difficulty  as  that  ot  NO 
due  to  their  coupling  through  the  NO,  cutalvtic 
cycle  Although  a onc-dimcnsional  model  can 
approximate  the  low  NO  column  measure 
ment  above  M>"N  In  using  the  propel  seasonal 
solar  flux,  temperature,  and  o/onc  (Noxott 
ct  ill.,  ll>7Saf.  a full  explanation  involves  the 
coupling  ol  three-dimensional  transport  and 
photochenustrv  tNoxon.  I *> 7 S > Nevertheless, 
on  the  average,  simple  one  dimensional 
models  do  produce  NO  and  NO  distributions 
that  are  consistent  with  stratospheric  data 


Stratospheric  NOi  abundance  as  esiima 
ted  In  photochemical  models  tails  tai  shoit  ol 
the  observed  amount  (Noxoii«7<i/  , l‘>',sbf  I Ins 
could  be  due  to  errors  in  the  accepted  kinetics 
rate  coefficients  as  suggested  In  Noxon  «7  ill 
tl*)7Xa)  But  this  involves  extiapolalion  ot  pie 
hminaiv  data  at  higher  temperatures  than  are  m 
the  stratosphere  f ortunatelv.  this  unsettled 
state  ol  attaus  is  ot  minimal  importance  to  this 
stlldv 

I he  current  estimates  on  stratospheric 
VO-  concentrations  are  all  consistent  with  the 
model  I he  upper  limit  ol  0. 1 ppbv  at  IS  ,1  km  In 
Murcras  «7  <;/.  (1978)  and  the  ' ppbv  at  '0  km 
deduced  bv  I vans  cr  ul.  (1978)  compare  well 
with  the  results  from  diurnal  calculations  utulei 
comparable  conditions  ol  0 0b  ppbv  and  2 1 
ppbv.  respectively 
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O('P): 

I he  local  coneentrations  ol  l)(‘P)  are  directly 
proportional  to  the  local  o/onc  concentrations 
and  solar  flux  intensity  . Any  disagreement  be- 
tween theory  and  experimental  data  on  O(V) 
can  be  directly  traced  to  disagreement  in  Ot. 
which  is  more  dependent  upon  atmospheric 
motions. 

a,  110  and  CIONOj : 

I he  onlv  av  ailahlc  data  on  l' I and  I'll)  are  I tom 
Anderson  <7  al  ( 1977,  1978)  As  the  discussions 
on  I'll),  indicated,  the  high  v allies  ot  I'll)  aie  not 
reconcilable  with  current  theoretical  undei 
standing  and  must  be  considered  the  ina.ioi 
weakness  in  current  theory  I he  low  I'll)  meas 
urements  are  in  good  agreement  w ith  theoiettcal 
models  I he  only  possible  exception  is  the  late 
ot  decrease  in  I'll)  as  one  appi ouches  the  nopo 
pause  One  dimensional  models  vield  a much 
more  gradual  decrease  I he  two -dimensional 
calculation  bv  Pyle  ( 1978)  shows  that  neai  W N 
(where  Anderson  <7  ill  (|977)  collected  then 
data)  there  indeed  is  a much  sliaipei  veitical 
giudicnt  trom  25  to  AA  km  I his  is  appaientlv  a 
direct  result  ol  the  detailed  liunspoit  piocesses 
in  this  region  It  should  be  pointed  out  that 
these  tw o dimensional  model  lesults.  although 
picliiiiin.il  v . aie  luglili  encouraging  \lthoiiph 
I'lONO-  was  shown  to  be  a cntical  species  m 
the  assessment  »»t  ll  M impact  on  stratosphciic 
o/otte.  up  to  now  onlv  an  uppci  limit  on  stia 
tosphene  abuinlance  has  been  suggested  I lie 
limit  ol  0 t ppln  at  20  k in  is  -it  least  a lactot  ot  2 
higliei  than  the  model  indicates  Hence,  it  is  an 
inconclusive  icsiilt.  although  somewhat  icns 
sin  mg 

HO.  IIOj  and  IfjOj 

I he  concentiations  ol  HO,  shown  in  I ig  -f  are 
all  diuinallv  aveiaged  values  Both  lit)  and 
lit'  undei  go  giadual  changes  tliiouglnuit  the 
dav  light  liv'd t s I oi  both  III)  and  III)-  the  lalio 
ol  noon  values  to  the  dnitnullv  aveiaged  values 
could  be  as  much  as  s m the  vvmtei  henuspheie 
high  latitudes  In  the  piesent  ease  ( <1*  tall 
ei|iunox)  the  ratio  is  appioxunatelv  2 A in  the 
altitude  lange  ol  10  20  km.  and  then  giaduallv 
dccieases  to  around  2 at  4s  km  Because  ol  tins 
large  averaging  tactoi  and  the  gradual  diurnal 
changes,  onlv  diutnal  models  should  be  used  to 
compute  vv  it  h expei  internal  data  Between  10  -It' 
km.  the  model-estimated  III)  piotile  at  com 
paiable  times  ol  the  vf av  and  tall  conditions 
tails  on  the  low  side  ol  the  data  bv  Anderson 
cr  ill  (l')’S).  the  computed  avetage  conoentin 
lion  being  t b \ It)  cm  ' -is  compared  to  an 
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avoiage  measured  value  ot  S \ It)  cm  1 At  '0 
km  the  local  \alui*  neat  sunset  is  about  a lactot 
ol  less  than  the  data  h\  Mihelete  ri  ill  t ll>”) 

I he  HO  dim  nal  pioliles  eompaie  \ei \ vv ell 
w ith  the  as  ad  it'le  data  between  '(>  40  kin  \t  the 
cotiesponding  tune  the  theoietieal  pi oliles  both 
(solai  zenith  angles  ol  41  and  SO)  tail  in  the 
middle  ot  the  data  tauge  (I  ig  <>'  1 Ins  is  pvob 
ablv  tin tuitous  sinee  there  is  no  eoi loborating 
data  on  lU,  II  t>,  and  o/one  which  are  essen- 
tial m determining  the  loeal  eoneentiations  ot 
lit'  rnlotiunatelv . in  the  entieal  legion  ol 

10  <0  km  there  aie  no  data  on  HO  and  lit' 

I he  onlv  estimated  nppei  limit  value  toi 

11  I'  at  km  is  well  above  the  model  eon 
initiations  to  be  ol  teal  use  I 'ue  to  the  eenttal 
i ole  plaved  hi  HO,  in  deteimmmg  the  ettcet 
ol  changes  in  NO,  eoneentiations  on  o/one, 
ineasuiemenls  ol  lit'.  Ill'  . and  II  O in  the 
lovvei  stratosphere  have  heeotue  more  nil 
poll, nit  than  evei 

Overall,  the  data  on  III),  seem  to  suppoit 
the  adei|iiaev  ol  the  tlieoiv  lot  the  piesent  appli 
eation.  but  lot  othei  piognostie  applications 
lelinement  on  the  lit',  budget  and  ehemieal 
ev eles  must  be  eat  tied  out 


2,2  Possible  \ ariations  in  Stratospheric 
O/one  Kelated  to  the  ll-\eur  Solar 
Cycle 

I he  suggestion  ol  a relationship  between  o/one 
and  the  solai  cvclo  dates  to  at  least  I 1 0 ( Humph  revs. 
I *4 10'.  and  Nations  reseatehets  have  tovmd  cot  vela- 
lions  ol  sunspot  numbei  with  total  o/one  (Willet, 
I'RO,  I’aet/old  a al , I47J;  (.’hnstie,  I**7-'.  Aiigollatid 
Korshover,  |4"’(')  Researchets  studviug  total  o/one 
lluetuations  lot  the  period  from  ll'5S  to  the  present 
have  in  general  found  less  evidence  to  support  an 
o/one  solai  correlation  (Angcll  and  Korshover, 
I'f’Na;  I ondott  and  Oilmans,  I')’1),  but  the  period 
over  which  o/one  was  studied  max  not  have  been 
long  enough  to  properlv  establish  a long-term  cycle 
Several  plausible  mechanisms  linking  the 
hvpothesi/ed  o/one  variations  with  the  solai  cycle 
have  been  suggested  Ruderman  and  Chamberlain 
( I *4 TS>  first  proposed  that  solai  modulation  ol  the 
source  ol  nitric  oxide  piodiiced  by  galactic  cosmic 
iavs  might  explain  the  o/one  record,  since  the  cosmic 
i av  source  ol  NO  is  approximatelv  h'1 , smallet  at 
solar  maximum  than  at  solar  minimum  A 
subscipient  studv  bx  Crut/enflRTMIailed  to  support 
this  mechanism  (see  also  Rudetman  *•/  .// . |4>) 
/eivtos  and  Crut/en  ( I ‘4 "■  5 > suggested  that  solai 


proton  events  max  also  prov  ale  a link  between  o/one 
variations  and  the  solai  cycle,  but  the  studv  was  not 
conclusive  Since  solai  pioton  events  aie  localized 
and  ol  short  duration  (Heath  < t ,tl . 1477),  it  is  not 
cleat  how  thex  max  ailed  global  o/one 
concentrations 

I oug  term  vauations  in  the  solai  uv  this 
pi  oxide  a thud  mechanism  ol  potentially  much 
gieatei  significance  lot  o/one  vat labihtv  in  a global 
sense  \V  e have  onlv  a few  sets  ol  measmements  that 
suggest  theie  max  be  a signilicaut  vanation  in  the 
solai  llux  between  ISO  and  140  nm  tiom  solai  maxi 
mum  to  solai  minimum  (Heath  and  Ihekaekaia. 
1477)  In  view  ol  the  possible  impot lance  ol  solai  uv 
variability  on  o/one.  we  shall  attempt  toestimate the 
elleets  such  a variability  would  have  had  on  the 
o/one  tccoiil  to  date  We  ate  pamculailv  coneeined 
with  the  o/one  record  neat  40  km  whete  both  the 
proposed  v a t labilitx  and  chlotme  increases  due  to 
past  tluorocarbon  telease  should  have  had  then 
greatest  elleets  I'allis  and  \ealv  (|4’S)  Inst  studied 
this  pioblem  with  a steady  state  model  without 
chlorine  chemist ix  Out  tune  dependent  model 
1 1 ut hot  1 1 ill . 147’’,  \ AS,  147b)  includes  a lull  set  ol 
cltloi  me  chemistry  tenet  ions  and  calculates  lent  pet  a 
lute  sell  consistewtlx  \s  ssc  shall  see,  out  calculated 
o/one  and  tempeiattue  responses  dillei  sigmticautlv 
Irom  those  calculated  bv  I'allis  and  \ealv 

I he  initiative  tiaustei  model  tot  the  tempeia 
tine  ealculation  used  is  described  in  I uihoi  <■/  til 
1 14 ’ 'I  I he  lull  set  ol  ehemieal  leacttons  is  given  m 
\ppendix  \ (lable  \ I)  Photolysis  tales  ate  com 
piitod  at  each  altitude  and  tune  step  consistent  with 
all  minot  species  disttihutions  and  the  specilied  solai 
conditions  these  calculations  note  made  using  the 
v eision  ol  t he  model  m vv  Inch  the  watei  v apoi  pi otile 
is  held  constant  l'onsei|iiently.  om  total  Hi',  at 
solai  maximum  and  minimum  is  somewhat  oxot 
and  under-estimated,  tespectively  I lie  elleets  should 
not  seriously  altet  om  main  results,  since  the 
calculated  local  vanation  loi  lit',  in  om  model  is  at 
most  ?()«', 

\ Mriahilily  in  the  Solar  l lux  and  Its  Kffect  on  O/one 

I igure  7 shows  the  ratio  ol  the  ultraviolet  (uv) 
llux  at  solai  minimum  to  the  uv  flux  at  solai  maxi 
ilium  as  a lunelion  ol  wavelength  (Heath  and 
Ihekaekaia.  I1)’  \ N ASA.  14”)  I he  measmements 
shown  bv  A and  represent  broad  band  photo 
mettic  observations  I hose  shown  bv  • vvete 
obtained  using  identical  double  monoeluomatois 
I he  measured  vaiiabihtv  below  7NN  nm  depends  on 
onlv  three  photoinetiie  obsei  rations  taken  neai  I St' 
nm  in  |4(vh,  |4(a',  and  l"’('  | he  llux  vanation 


I K.  7 Ratio  o(  ullmv  inlet  (lux  at  solar  minimum  lo  the  flux 
al  solar  maximum 

measured  with  these  instruments  at  longer  wave- 
lengths overlaps  the  region  measured  hv  the  double 
monochromators,  the  photometers  show  about  a 
IV’,  larger  variation  than  the  double  monochroma- 
tors We  believe  the  double  monoehroniatot  meas- 
urements at  these  wavelengths  are  more  neatly 
correct,  because  they  were  marie  with  identical 
instruments.  Since  there  are  no  monochromator  data 
below  255  tun.  we  applied  a correction  ol  about  I 5' , 
to  the  average  ol  the  photometric  measured  flux 
variation  neat  ISO  tun.  I he  curve  in  fig.  7 marked 
"II"  thus  represents  what  we  consider  to  be  the  largest 
probable  variation.  I lie  curve  marked  “A"  gives 
more  weight  to  the  measured  variation  ncai  755  inn 
and  is  also  consistent  with  balloon  measurements 
which  indicate  that  the  lluv  variation  between  700 
and  770  ntn  was  not  larger  than  10*',  trom  1977  to 
|9"(v  (Simon.  1978).  Hie  curve  marked  "( repre- 
sents an  intermediate  case  and  is  used  lor  most  ol  the 
calculations  shown  in  the  following  sections 

I he  calculated  o/one  response  in  our  model 
depends  fairly  critieallv  on  the  relative  flux  variation 
.dune  and  below  755  tun  I’hotons  with  wavelengths 
less  than  755  tun  can  add  to  the  total  atmospheric 
reset  s on  ol  odd  oxygen  according  to  the  i eaction.  O 
r hr  -(HO  I’hotons  in  the  entire  range  ol  wave 
lengths  considered  here  are  able  to  slult  the  balance 
between  at  tunic  ox  v gen  and  o/one  according  to  O,  t 
lit  • O t (>  . but  bv  lai  the  most  efficient  are  those 
with  wavelengths  longer  than  755  inn  I he  point  is 
illustrated  m I ig  N I he  curves  markcil  VI  and  Ml 
show  the  calculated  o/one  increase  when  the  solai 
tlux  is  only  varied  below  755  ntn  by  the  amounts 
indicated  in  curves  V and  M ol  I ig  7.  respect ivelv  V 
larger  llux  at  these  wavelengths  produces  more 
o/one  I he  curves  marked  A 7 and  H7  show  the  o/one 
response  to  llux  variations  above  755  tun  onlv  V 
large  llux  at  these  longer  wavelengths  slults  the 
balance  of  odd  oxygen  m lav  or  ol  atomic  oxygen  so 
that  o/one  is  decreased  I he  net  result  lot  llux  vana 
lions  ov ei  the  entire  wavelength  lattge  is  to  increase 


l'l<*.  S.  Permit  change  in  (I  from  mitt  miniimmi  lo  xotar 
maximum  |(niax-niin)  min  s 1 00 1 the  Irltrrx  V amt  H refer  lo 
the  dux  change  runes  giien  in  l ig.  I.  l or  rases  VI  ami  HI  the 
solar  flux  wax  onlv  \arietl  hrlou  IW  on,,  l-or  raxes  \2  ami  H2 
the  solar  nux  was  onlv  varievl  above  2SS  mo 

o/one  ovet  most  altitudes,  as  show  n bv  the  curves 
marked  M and  A.  and  the  total  o/one  response  is 
approximately  the  sum  of  the  variations  calculated 
separately  I Ins  demonstrates  the  nee.)  lot  a precise 
determination  ot  the  relative  flux  variations  in  these 
wavelength  bands.  Since  the  net  o/one  response 
depends  on  the  superposition  ol  two  eflects  ol 
opposite  sign,  it  is  cleat  that  the  altitude  ol  the 
maximum  response  depends  on  this  relative  llux 
variation. 

Kffeetx  of  t emperature  l eedhack 

1 he  primarv  source  ol  heating  toi  the  upper 
stratosphere  is  the  absoi  pi  ion  "I  solai  photons  in  the 
1 1 art  lev  (700  100  mill  and  Muggins  ( 1()()  nm)  bands 

ol  o/one.  We  use  the  method  ol  I ind/en  and  W ill 
( 197.1)  to  calculate  the  heating  rate  liom  these  bands 

I igure  .1  shows  the  calculated  change  hi  temperature 
and  o/one  when  the  solai  llux  is  perl uihcd  from 
minimum  conditions  bv  the  amount  shown  bv  curve 

II  in  fig  7 Out  calculations  show  a temperature 
inciease  ot  no  nunc  than  about  I ' degrees  tieai  45 
kin  I he  I iinl/eii  .nixl  \\  ill  paiametei i/ation  igmues 
llux  perturbations  below  7.''  5 nm.  but  since  these 
wavelengths  contiibutc  no  mote  than  about  10* , ol 
the  ambient  heating  i.itc  at  4s  km  out  calculated 
’em pei  atm  c change  shou  lit  be  coi  i ect  tx»  w it  Inn  1(1* , 

I he  lelt  panel  ol  I ig  ' shows  the  ellccts  ol  including 
tenipeiaiuie  feedback  on  calculated  ti  conceutia 
lions  \\  it li  tenipeiaiuie  feedback  the  |'cak  i espouse 
ol  o/one  is  reduced  due  to  the  tenipeiaiuie  depend 
dice  of  the  reactions  t)  I O,  •()  t O andO  * O I 
VI  • (),  t M.  and  the  altitude  above  which  o/one 
decreases  is  lovveicd  liom  50  km  to  41  km 

Om  calculated  tempeiature  increase  at  solai 
maximum  is  significantly  smallci  than  that  calcu 
lated  bv  I'allis  and  \ealv  (|9'S)  I Ins  lemains  true 


even  lot  assumes)  tlu\  poiim  bations  as  huge  ;»n  then 
case  I We  believe  then  pitman  t henna  I icsponse  is 
vine  tv'  heating  b\  wavelengths  beh'vv  .'It'  mil.  Mine 
an  examination  ol  then  ealenlatevl  tenipeiatine  te 
spouse  tin  ease  I tsvinbi'l  v>t  then  papei'  slnivvs  a 
Mgmtieant  meiease  in  tenipeiatine  (.'OK  neai  4'  km' 
even  though  the  xolai  this  is  not  pert  in  I'evt  ahov  e .'It' 
inn  anti  then  ealenlatevl  ozone  eotieeniiation  at  that 
altitnvle  is  vleeieasevl  (eoinpate  I igs  .'  ami  ' ol  then 
papei'  \ vletailevl  ealeulation  vil  the  ambient  heating 
rate  ttoni  ('i  absv'tptivin  at  4s  km  shows  that  the 
heating  hv  wavelengths  below  .'It'  nm  is  onlv  0 O'1 , 
ot  that  lot  the  entile  interval  liom  I ’s  to  'M'  nm,  sv' 
that  this  pei tin hations  in  tins  legion  shvnilil  not 
gieatlv  atleet  stiatv'spheue  tenipeiatine  I ignie  *> 
shows  that  tenipeiatine  leeilbaek  is  impintant  tv'i 
ealenlating  the  eliange  in  o/otte  above  abviiii  km. 
but  the  ealenlatevl  tenipeiatine  eliange  is  not  ne.nlv 
so  sensitive  tv'  tins  vaiiatn'iis  as  that  ealenlatevl  hv 
C'a Ills  ami  Nealv 


Comparison  With  the  Ozone  amt  temperature 
Reeortls 

I he  ealeiilatnvns  vleseubevl  thus  lai  vveie  pei 
lormevl  using  a steavlv  state  vet  sum  ol  >'in  imulel 
lleie  we  use  om  tune  vlepemlent  movlel  tv'  examine 
the  atniv'sphei n tvsponse  tv'  a time  vaiving  solai 
(lux  W e iisevi  tlw  release  rates  lot  t I t I.  ami  I I i I 
v't  MeC'aithv  er  a.  1 1 4 "'  ami  the  solat  this  v ai  latn'ii 
l' in  I ig  ' to  ealeulate  a tune  Instv'iv  loi  ozone  ovei 
the  past  4'  veais  I he  this  vaiiatn'ii  li'llv'vveil  a 
sunisoivt.il  II  veai  evele  with  maximum  ami  mini 
mum  eon  espousing  tv’  the  veais  l'"('.  I'M".  |4ss, 
|4n4  ami  |44.'.  I*"!,  |')(>.|.  |0's  iv-speetivelv  I he 
ealeulatu'ii  melmlevl  all  v't  the  leevlb.u  k ineehanisms 
vlisv  tissevl  abv’v  e ami  simulates  avei  age  eoml  it  ions  toi 
'('  latitmle  Seasonal  v anal  mils  in  the  solai  llux 
vveie  not  inehntevl.  smee  all  v'ne  vliiuensn'iial  niv'vlels 
tail  tv’  aeematelv  simulate  the  seasv'iial  behavn’i  v't 
v’/v'iie  vine  tv'  then  laek  v't  hon/ontal  ttanspvut  We 
eotnpate  vmii  ivsiilts  to  the  analvsis  bv  Vngell  ami 
Kv'islu'vei  (|4'S.i.b,e)  m vvliieh  seasv'iial  etteets  ate 
partlv  leiimveil  bv  plotting  seasv'iial  vtev  lativ'tis  Invin 
seasv'iial  means 

I he  ealenlatevl  pereentage  vleviatn'ii  v't  total 
v'/v'iie  1 1 om  the  mean  is  shown  in  I ig  It'  loi  the 
pet  iv'vl  I *>s'  to  1 4 'S  along  w till  ( he  I v't  a I ozone  leeoivl 
tv'i  II  stations  m Nv’itli  Vmeiiea  amt  14  stations  m 
I mope  t Vngell  ami  hoishover,  |4’Sa'  I he  iiVvitvl 
liv'iii  North  \metiea  ami  I mope  was  eliosen  tv'i 
ev'inpaiison  beeause  these  aieas  have  the  hugest 
niunboi  v't  siativ'iis(v'thei  than  the  Sv'viet  I niv'ii'aml 
aie  eonsuleu'vl  to  be  the  niv'st  .urinate  ami  tellable 
I he  ertv'i  bats  tv'i  the  v'bsetvativ'ns  ate  two  simulant 
vlev latu'iis  v't  the  mean  It  sluuilvl  he  pointed  v»ut  that. 


as  slu'wii  in  lag  4.  the  hugest  ehanges  in  evuieen 
latum  oeein  above  '.'  km  vvlieie  the  eheinieal  tune 
n'lisiaiit  tv'i  v'vlvl  v'xvgen  is  v'ulv  v't  theonlei  ol  weeks 
tv'  months  i v'lisevpientlv . the  lesiilts  in  I ig  It'  sliv'w 
v ei v little  tune  vlelav 

l leailv.  v'iii  ealenlatevl  v'.’vme  vaiiatnm  tv'i  this 
pei u'vl  is  at  least  inn  iiieonsistent  with  this  tceonl 
the  obseiveil  v'Zv'iie  meiease  timing  the  hue  |4(.ii's 
appeals  tv’  he  teasonahlv  evutsistent  with  out 
assiuncvl  eliange  in  solai  llux  We  ealeulate  an 
meiease  in  total  Oi  livun  the  solai  minimum  in  I4(<.j 
tv’  the  solat  maximum  m |4<.4  ot  x |> , | ,'i  sunilai 
v ai  tat  ii'ii  m t he  solai  llux  bet  ore  the  nm  otluet  ion  ol 
Cl  M's  we  I'htam  an  meiease  ol  ' V,  1 Ins  mav  he 
ev'inpaiVvi  t>'  the  obsetvevl  V,  meiease  in  totalOi  lv>i 
the  Noithei  u llemispheie  I tom  the  eat  Iv  (>t''s  i,’  1 4 'il 
ami  meiease  tv'i  tv'talt'i  tv'i  the  Southern  Ileum 
plieie  tiv'in  14(v|l  tv’  |4<<s  We  ealeulate  a vleetease 
tiv'in  1 4 *|I  tv'  1 4 * ' v't  .'  S',  w hit'll  agiees  well  w it  h the 
ot'seivevl  I , vleeiease  m the  Nv'ithetn  Menus 
pliete  I he  vleeiease  m ohseivevl  ozone  appeals  U’ 
Ii.iv  e lev  elvul  ot  I allot  1 4 but  we  ealeulate  a tin  thei 
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decrease  through  1475.  Our  calculated  decrease  is 
primarly  due  to  the  solar  llux  variation,  since  the 
1975  minimum  is  only  I ..V’;  smaller  than  that 
calculated  without  the  influence  of  increasing 
chlorolluoromethanes. 

Clearly  there  are  local  disagreements  between 
theory  and  observation  throughout  the  period 
show  n.  Fv  identic  some  processes  that  control  o/one 
on  shorter  time  scales  (3-4  years)  are  not  accounted 
for  by  the  model,  and  the  solar  u\  connection  and  the 
effects  of  increasing  chlorolluoromethanes  w ill  only 
be  revealed  by  looking  at  the  general' long-term 
picture.  Thus,  the  apparent  discrepancy  between 
theory  and  observation  after  1472  does  not  neces- 
sarily invalidate  the  background  o/one  cycle  caused 
by  solar  uv  llux  variations  and  chlorine  changes. 

I he  calculated  o/one  deviation  for  the  levels  32- 
46  km  is  compared  to  Umkehr  measurements ( Angell 
and  Korshover,  1978a)  in  big.  II.  Again,  the  two 
regions  with  the  largest  number  of  stations  were 
chosen  for  comparison.  Umkehr  measurements  are 
much  more  difficult  to  make  than  measurements  for 
total  o/one.  I he  observations  should  accordingly  be 
viewed  w ith  caution  (sec  Putsch  and  I ing.  1973)  At 
these  levels.  Angell  and  Korshover  estimate  there 
was  an  8r(  increase  in  north  temperate  latitudes 
between  1962  and  1973,  and  a slight  decrease  after 
1974  ( I lie  observed  decrease  may  be  fictitious  due  to 
the  presence  of  aerosols  in  the  stratosphere  after  the 
eruption  of  Fuego  in  1974.)  We  calculate  a 5.8r, 
increase  for  this  layer  from  the  solar  minimum  in 
|4(vj  to  the  maximum  in  1969,  and  an  S.2‘ , decrease 
from  1 969  to  1975.  I his  may  be  compared  to  an 


increase  of  6.8*';  Irani  solar  minimum  to  maximum 
and  a decrease  of  6.4%  from  maximum  to  minimum 
for  the  pre-CFM  atmosphere.  As  shown  in  Fig.  1 1. 
our  calculated  decrease  begins  earlier  and  is  larger 
than  the  observed  decrease.  I’his  difference  between 
the  model  and  observation  may  be  indicative  of  the 
presence  of  other  physical  processes  controlling  the 
local  o/one  concentrations  that  are  not  included  in 
this  model  (see  London  er  <r/..  1977). 

As  Fig.  12  shows,  our  calculated  decrease  after 
1969  is  consistent  w ith  the  o/onesondedata  for  lower 
altitudes.  I he  Umkehr  data  for  32-46  km  is  shown 
again  along  with  our  calculations  and  the  o/one- 
sonde  record  for  the  layers  24-32  ktn  and  16-24  km 
(Angell  and  Korshover.  1978c).  Our  calculations 
agree  fairly  well  with  the  record  for  F'urope  for  these 
lower  layers.  Angell  and  Korshover  estimate  a 4<  < 
decrease  in  each  of  the  two  lower  layers  as  an  aver- 
age for  the  north  temperate  latitudes  as  a w hole  We 
calculate  a decrease  front  the  solar  maximum  in  1969 
to  the  minimum  in  1975  ol  7.S' , for  the  24-32  km 
layer  and  4.8('i  for  the  16-24  km  layer. 

f or  the  layers  shown  in  f ig  12.  we  calculate 
decreases  of  o.4‘ , . 6 41 , . and  4.6‘V . respectively . for 
the  same  solar  llux  variations  front  the  maximum  in 
1969  to  the  minimum  in  1975  without  the  influence  of 
increasing  chlorolluoromethanes.  For  this  period 
changes  in  o/one  due  to  the  release  of  CFM's  are 
apparently  small  compared  to  changes  resulting 
front  our  assumed  solar  flux  v ariation,  rhereforeany 
effort  to  detect  a trend  in  o/one  due  to  the  ('I'M 
releases  must  take  into  account  the  effects  ol  solar  uv 
variability 
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t It),  lit  omparison  of  ralrulatrd  lavrr-mran  o/onr  variations 
to  Ihr  analysis  of  l mkrhr  ilata  hv  \n*rll  amt  korshovrr  ( I »7Ma ). 
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I It).  12.  ( omparison  of  ralrulatrd  lavrr-mran  o/onr  variations 
to  Ihr  analvsis  of  o/onr  data  hv  Vngrll  and  korshovrr  (IWIHr). 


figure  I J compares  out  calculated  tempeiature 
variations  li»i  the  laseis  2b-J5  Km.  Jb-45  Km.  and  4<> 
55  Km  to  rocKetsonde  temperature  data  t \ngcll  and 
korshoset.  Id'Sb)  l he  lempetatuie  maxima  and 
minima  appear  to  he  nearly  m phase  with  the  solat 
cycle  sanation  except  lot  the  highest  laser  helot e 
ldr>S  (although  the  equatorial  stations  show  .1 
decrease  tluough  the  entire  tune  interval.  see  \ugcll 
and  Korshoser.  Id?8b)  We  calculate  a decrease  in 
temperature  from  ld<W  to  Id's  ol  S 5 k.  o 5 k.  and 
I I k lot  layeis  40-55  Km.  Jo-45  Km.  and  15  Km. 
lespeclisely.  in  reasonahle  agreement  with  Mtgell 
and  korshover's  estimate  ol  about  a 5 k deciease  lot 
the  middle  and  upper  stratosphere  alter  Id’O 
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H(».  I.r.  t < mi  pah  si  mi  of  calculated  laict  mean  trm|icrlurr 
changes  lo  Ihr  analssis  of  rmltrlsiimlr  data  In  Vngcll  anil 
korshoif r (|*7M»). 


Nummary  and  Conclusions 

I he  o/one  response  to  a variety  ol  solat 
conditions  was  examined  and  shown  to  depend  on 
the  relative  ratio  ol  u\  sariabihts  above  and  helow 
about  255  tun  I he  temperature  increase  due  to 
increased  heating  b\  o/one  absorption  significant l\ 
affects  the  total  o/one  variation  above  U>  Km  O/one 
changes  associated  with  solat  this  sanations  ate 
prinvartlv  due  to  changes  in  the  phs'tolysis  Kites  lot 
O and  Oi.  and  only  second. mis  due  to  changes  in 
other  minoi  species  W e used  w hat  we  consider  to  be 
the  most  reasonahle  sanation  in  the  solar  us  flux 
consistent  with  available  data  to  calculate  a tune 
history  lor  o/one  and  temperature  We  conclude  that 


the  o/one  rccoid  as  published  bs  Vugclt  and 
koishovci  (Id'ga.c)  is  apparently  consistent  with 
solat  flux  s .illations  ot  the  magnitude  used,  although 
more  data  are  required  to  lulls  establish  this 
lelatiouship  Some  ol  the  less  satislactois  icsults 
abv'se  '2  Km  requite  luithei  analysis  and  moie  quan- 
titatively reliable  data  We  calculate  an  o/one 
deciease  ol  Jo<!t  lor  the  '.’-40  Km  lasei  liom  solai 
aseiage  conditions  in  |dt>n  to  asetage  conditions  in 
ld77  due  to  the  release  ol  t’lO,  liom  lluotocat bons 
In  our  view  this  decrease  would  be  ditticult  it  not 
impossible  to  detect  tront  the  given  o/one  iccoul  it 
solai  I In  x sanations  ol  the  magnitude  considered 
here  were  present  dui  mg  this  petiod  because  the  local 
o/one  sanations  caused  bs  solat  flux  sanations 
could  be  as  large  as  I J1',  Insolat  as  the  piesent  calcu- 
lations should  be  considered  tentative,  a long  tune 
senes  ol  measuiemcnts  ol  the  us  (lux  as  a (unction  ot 
wavelength  is  needed  to  further  establish  the  sahdits 
ot  these  cttects 


2.3  Kffccls  of  a Solar  Kclipst'  on 
Stratospheric  ( licmist rv 

Solai  eclipses  are  Known  to  allect  measiiiable 
quantities  in  the  uppei  atmosphere  Various  studies 
have  discussed  the  response  of  ionospheric  election 
and  positive  ion  densities  to  an  eclipse  (I  atidm.it K «•/ 
.;/ . Id'O.  Marriott «-/ «//. . Id'7,  -\nastassiades.  Id’O) 
Other  studies  have  discussed  the  possibility  ol 
changes  in  stratospheric  and  mesosphenc  o/one 
during  an  eclipse  Sti.ui/  (Idoll.  using  a Hobson 
o/onometei.  measured  appioxunatels  a4\  met  ease 
in  total  o/one  shortly  allot  the  maximum  phase  of  an 
eclipse  m which  only  about  SO1,  ol  totahts  was 
reached  Hunt  ( |d(i.si  attempted  to  explain  this  total 
o/one  met  case  bs  using  I'hapnian  reactions  in  a 
theoietic.il  studs  ol  the  effect  on  o/one  resulting 
from  an  eclipse  Hunt  found  that,  at  most,  an 
increase  ot  onls  0 b‘\  is  to  be  expected,  and  that 
o/one  should  be  affected  only  above  4'  Km 

In  latei  sibsersatums.  Kandhassa  (Idost  and 
Haltaid  <■/  ill  ( Idpdy,  both  using  locKet-boine  o/one- 
sondcs.  lound  measurable  incivases  in  o/one  densitv 
above  50  Km  lot  the  Nosembei  Idoo  eclipse  neat 
S Kandhassa.  tot  example,  measured  a tactoi  ol 
Jo  increase  in  o/onc  at  5'  Km  relative  to  the  das 
beioie  Howes ei . Kandhassa  (Id'J)  found  no  signiti- 
cant  o/one  change  between  50  and  55  Km  lot  the  July 
ld'2  eclipse  at  t>5  \ tiiasiucK  <•;  ill.  ( |d’4(  lound  no 
change  m total  o/one  using  a Hobson  o/onomctei 
tin  the  Match  Id'O  eclipse  I lies  condudeil  that 
pies  ions  measureinents  had  not  collected  loi  solat 
limb  dai Kemng  However,  Oshcrovich  <•/ .;/  (|d'4t. 


lor  the  luls  I'f’’  eclipse,  claim  to  ha\e accounted  lot 
limb  darkening  and  set  lind  apptoxtmatels  a V\ 
increase  in  the  total  o/onc  column 

l leai  Is . inconsistencies  exist  in  the  obxci  s at  ions 
and  in  then  utilization  lot  s ululating  ozone  lelated 
atmospheric  chemistis  as  presently  undeistood  In 
addition  to  ozone  ohsersations.  measurements  ol 
othei  nunoi  cs<nstituents  duiing  a solai  eclipse  could 
ptosule  salutation  ol  the  short  lifetime  chemistis 
used  m atmospheric  models  Consequentls . experi- 
ments loi  upcoming  solai  eclipses,  when  pioperls 
supported  hs  theoietical  analysis,  could  contribute 
significantly  to  piesent  undeistandmg  ot  atmos 
phene  chemistry  In  tact,  propei  measuiements  ot 
trace  species  concentrations  during  an  eclipse  could 
proside  a direct  demonstration  that  curtentls 
proposed  NO.,  HO>.  and  I'lO.  catalsttc  cycles  ate 
indeed  concurientls  tunctiomng  m the  stiatosphere 
in  the  manner  suggested  bs  laboiatois  chemistis 
While  detailed  measurements  ot  diurnal  sanations 
could  piosule  sunilai  knossledge.  the  ditteience  m 
tune  scale  betsseen  the  diurnal  cycle  and  a solai 
eclipse  esent  suggests  tliat  the  lattet  esent  mas  pio 
side  a clearei  pictuie  lot  understanding,  as  ssill  be 
discussed 

I he  put  pose  ot  this  souls  is  to  examine  ilteo 
reticalls  the  expected  etlect  ol  a solai  eclipse  on 
stratospheiic  minor  constituents  I’nmais  emphasis 
is  gisen  to  calculations  m ads  a nee  ot  the  total  eclipse 
tliat  ssill  occur  osei  North  Xmeitca  on  t ebmais  .'o. 

I *•'’**  (t-'iala  and  I ukac.  I 7 "■  > Variations  sunilai  lo 
those  ci'iuputed  tot  this  particulai  case  slumld  be 
expected  tor  othei  total  eclipses 

Methodology 

I he  III  one  dimensional  model  ol  tiopos 
phene  and  stratospheric  transport-kinetic  piocesses 
liom  0-5.'  km  ssas  used  in  these  calculations  I he  lull 
set  id  chemical  reactions  is  gisen  in  \ppendix  A 

I I able  A lt 

l ifts  degrees  north  is  the  latitude  at  which 
totality  ssill  be  the  longest  (5  minutes)  lot  the 
I ebruars  I '>"’**  eclipse  In  attempting  to  simulate 
tspical  expected  atmospheric  conditions  tot  the 
I ebruars  'b.  I s» ’sJ  Sl>|.,|  eclipse,  sse  iiuulitied  the 
model  bs  setting  the  temperature  ft  ouis,  I'Od). 
tropospheric  water  sapoi  (Oort  and  Rasmussen, 
I ‘Tit.  and  the  ozone  distribution  (Wilcox  and 
Belmont.  I*)7?)  at  50  N w inter  conditions  I atomic 
and  the  solar  declination  angle  in  the  model,  neces 
sats  to  calculate  the  solai  zenith  angle  lot  photodis 
socialism,  were  set  to  50  and  -S  S '.  respect  is  els 

While  ozone  in  the  model  was  held  tixcd  lot 
most  calculations  in  ordci  to  simulate  lebtuats 
50  N conditions,  the  model  ssas  also  i tin  ssith  ealeu 
lated  ozone  to  examine  the  expected  tesponse  ot 


ozone  to  an  eclipse  Oui  analssis  lias  shown  that 
fixing  the  ozone  distribution  dm-s  not  siguificantls 
allccl  the  tempoial  sanations  calculated  toi  othei 
species  during  the  eclipse 

Assuming  fixed  total  odd  nitiogen  and  odd 
chlorine  as  calculated  lot  mid  latitudes,  the  model 
ss  as  i tin  dun  nails  to  cs|Uilthi  nun  I lie  calculated 
daytime  total  column  lot  NO-  at  50'  N ol  I *)  \ Id1' 
molecules  cm  is  consistent  ssith  the  obseisations  ol 
Noxon  ( I sJ 7S) 

Solat  llux  sanations  dining  the  eclipse  sseie 
based  on  Hunt  (|4b5)  I he  eclipse  ssas  assumed  to 
stall  at  4 10  a m and  end  at  I'  IK)  noon,  svitespond- 
ing  to  the  I ebi  iiat  s I1!  ’*>  eclipse  toi  50  N m southei  u 
tannd.i  lotalits  ss as  assumed  tooccui  lot  I minutes 
beginning  at  10  4'  '0  a in  Based  on  measuiements 
ol  sks  bnghtness  dining  an  eclipse  (\ el.isi|uez.  |4’| 
and  I >aiidek.ii  and  1 ui tie.  |47|),  the  solai  flux  at 
lotalits,  ssas  i educed  to  10  * ol  the  unobscuicd  llux 

Results  and  Discussion 

I hose  species  basing  chemical  lilciimes  less 
than  a less  houis  should  be  expected  to  sais 
xigmttcuiitls  Itotn  nonnal  diunial  behasioi  during. i 
si'lai  eclipse  In  this  studs,  sse  ssill  emphasize  the 
sanations  expected  toi  those  species  most  mipottant 
to  the  chemistis  «'l  stiatosphenc  ozone  1 ocal  con 
cent  i at  unis  ot  the  species  could  be  ijiute  sai  table,  and 
theieloie  sse  should  locus  on  i el.it is c ellects  in  the 
calculations  i.ithei  than  then  absolute  magnitudes 

I he  model  calculated  i espouse  ot  ozone  dm  mg 
an  eclipse  essentialls  agiees  ssith  Hunt  ( |4<i5l  \ sig 
mticant  metease  in  O.  is  to  be  expected  m the  uppei 
stiatospheie  and  in  the  mesospheie  due  to  the  con 
set  sum  ol  Of  IM  u>  ozone  tluough  the  reaction  ()(  I’ I 
t O'  t M • if  r M accompanied  bs  decieased 
photolssis  ot  if  and  If  I lie  maximum  tnciease  m 
If,  lound  at  the  end  ot  totality , ssas  computed  to  be 
14  and  45> , at  50  and  NS  km,  icspectisels  1 arget 
peicent  changes  should  be  expected  in  the  meso- 
sphere Since  most  ol  the  atmospheric  ozone  is  at 
losset  altitudes  in  the  stiatosphere,  an  insignificant 
change  in  the  total  ozone  column  is  to  be  expected 

f igure  14  shosss  the  sanations  expected  dut mg 
the  eclipse  in  the  concentration  ol  nituc  oxide.  NO, 
and  nitiogen  dioxide,  NO-,  at  altitudes  ot .'(',  '(',  and 
40  km  I he  solid  line  indicates  notmul  diurnal 
behasioi  ftom  S a in  to  p nt  , while  the  dotted  line 
shosss  the  change  in  concentiation  expected  during 
the  eclipse  l he  decieased  solai  tlux  results  in  the 
NO  phi'tolssis  rate  decreasing  during  the  eclipse 
Also,  the  rapid  conseision  ot  O('l’)  to  if  dect eases 
the  importance  ol  NO-  * 0(1’)  -NO  * O-  I he  NO 
is  then  lapulls  conseued  to  NO-  dining  the  eclipse 
pnmaiils  bs  the  leactions  NO  > if  • NO  * If  and 
NO  * Ollt  • NO-  t ( I Bs  the  end  ot  totahts,  osei  a 
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f It;.  14.  ( alv'uUtrd  mpcinm  of  NO  and  Nt>  In  a ,ol«r  nli|iM-  Ihc  solid  linr  indicate  normal  diurnal  brharior.  I hr  dotlrd  linos  ,ho>*s 
(hr  r\prr(rd  vilar  rrlipsf  effect . 


factor  ot  If)  decrease  is  calculated  lot  NOal  20  and  >0 
km.  with  a factor  ol  tv 5 decrease  at  40  km  I he 
maximum  NO  concentration  also  occurs  at  the  end 
ot  totalitv  with  rallies  ranging  from  a laetot  ot  I 4 
higher  than  normal  at  20  km  to  h 4 at  40  km 

1 \cept  in  the  upper  stratosphere  (i  e . 40  km), 
the  NO  increase  is  not  direeth  proportional  to  the 
NO  decrease  1 his  is  due  to  the  relativeh  slow,  vet 
significant,  tormation  rates  toi  1'IONO-  and  VO. 
from  NO  I he  peak  increase  in  thel'IONO  concen- 
tration occurs  approximate!)  15  minutes  after 


totalits  due  to  the  slow  tot  illation  rate  llTO  * NO-t 
l IONO  is  increased  h\  20' , at  20  km  and  42'',  at  40 
km  NO.  concentration  increases  b\  ll‘  , at  20  km 
and  20' , at  40km  lheslow  reconversion  ot  I'lONO- 
and  N O.  to  NO  causes  the  NO-  to  tall  below  the 
normal  diurnal  concentration  at  the  end  ol  the 
eclipse  I he  NO  to  NO-  conversion  is  not  aftcctcd  in 
the  uppei  stratosphere  where  neithei  I'lONO-  not 
N O,  is  important  as  a temporal')  sink  toi  NO 

I he  expected  variations  toi  t'l  and  I'll)  are 
shown  in  I ig  15  Rapid  conversion  ot  t'l  to  l it) 


P 


on i mg  (he  eelipie  oeeun  pritnaiili  thtough  l I * U 
•CIO  * O I he  net  etti-el  on  the  CIO  eoneenitaiii'ii 
n .1  ileeieaie.  huttcu’i.  Jiic  1 1*  iiii'ii- 1 .1  pul  li>i m.ilii'ii 
ol  CIONO  n-Miltiilg  1 1 i'm  uieieaieil  NO  .nut 
rciluccil  CIONO  photoli io  I hi-  iimnmum  in  I'll' 
oeeun  appioiiniateli  '■>  inuiulei  altei  totaliti  I ln\ 
linn-  lajt  n eatm-il  h\  i ho  lion  It'i ni.it ii'ii  i .no  ol 
CIONO  \l  I 111-  M'l.ll  lhl\  llll  ll-.IM-N  .llll-l  ll'l.lllll  . I In- 
i'll' ei'iieenttation  iiteiean-i  line  in  met  eased  pln'to 
li  m\  ol  CIONO  1 ntle  eltaniii-  in  i'll'  n oipeeteil  .it 
40  kin  iilieie  CIONO  tinniiilton  n not  inipoit.int 
Noth  i'll  .mil  Ill'  ihonlil  decrease  in  I'l'iieen 
1 1. n ion  ilnnni;  an  eclipse.  .i'  ilioun  in  I i>:  !t<  i'll 
i-oiiei-iitiatioiii  nete  u-ilneeil  lelatue  to  the  itoini.il 
iluiiii.il  eoneentiatiotii  hi  a laetoi  ? I'  0.  anil  > *> 
ai  '0,  '0.  ami  40  km.  leipeeineli.  null  a ininiinnin 
oii-niiini;  at  the  eiul  ot  totaliti  Smallei  ileeieaiei  ate 
ealenlateil  toi  HO  . null  the  uiinitmitu  in  III' 
oeeui i unt  uppioviniutelt  Is  nnnutei  attei  totaliti 
With  the  ileei eau-il  unpoitanee  ot  l't  I’t  ami  NO 
ilnnni:  the  eelipie.  a halanee  hetneen  OH  ami  III'-  n 
inaintameil  hi  the  tapnlli  oi'enri mg  reaction  OH  * 
O • III'  * O aml  the  iluthtli  lionet  ii-aetion  lit' 

* Oi  -I'll  > .'O  Dili)  hi  ill  ojien  (OH  ami  HO  I Ion 
oeeui i at  a ilonei  late  tluoui;h  eoiiieision  to  iiatet 
(I'll  • HO  I.  H O (HO  * III'  'ami  UNO  (i'll  • 
Ni'  t \l'  Insijtiiitieant  ehatiiiei  in  the  eoneentia 
tiom  ot  H O ami  IINOi  neieealeulateil  \t  40  km,  a 
ihjtht  ineiean-  in  III'  n expeetfil  at  totaliti  line  to 
the  laiite  t a put  eoiiietsion  ol  I'll  to  IIO‘  loiloneil  hi 
the  ilonei  III'  Ion  iiieehaiinmi 


SumniHri  ami  t omluiioin 

I he  i emit i liom  thii  itmlt  nifijieit  that  ngniti 
i .i nt  ami  ileteetahle  lanatioin  .lie  eipeeteil  lot  lome 
ol  the  uupoi taut  itialoiphene  mtnoi  eonititueiili 
ilutiti)!  a n'lai  eelipie  Nueli  ohieiiatioiu.  p.utie 
ulaili  iimultaiieoui  ohienations  ol  tiaee  ipceiei. 
noulil  ileinoiiitiate  eleaili  the  Miiiultatu-ons  tune 
tioning  ol  the  latioui  iinpoitant  photoehi'imeal 
eatalitie  eielei  m the  stiatoipheie 

In  pi  meiple.  a umilai  ti  pe  ol  mtoi  illation  ean  he 
ol't.iiin-il  hi  ohieiiatioiii  iluiinj;  a iiomial  ilmin.il 
ei  ele.  eipeeialli  Oui  mi;  nun  ne  ami  mmet  Ilonei ei , 
ohieiiatioiu  ilut ui(:  a lolai  eelipie  otlei  n-ieial 
.nil antaiiei  oiei  that  ol  the  iliutnal  eiele  1 mi.  ai 
nai  iliieuni-il.  the  telatneli  ilioit  tune  ilutation  toi 
the  eient  loiloneil  a moie  ele.u  ulentilieation  ol 
n'lne  ol  the  maioi  teaetioin  in  the  uupoi  taut  eata 
line  evelei  1 he  lame  ilioit  time  iluiatii'ii  t?  > In' 
alio  muiuni/ei  the  intluenee  ot  atmoipheiie  tiam 
poit  pioeen  m altentiji  the  loeal  tiaee  ipeeu-i 
eoneentiatioui  \ lougei  ohieitation  penoil  noulil 
mhieet  the  loeal  eheinutii  to  the  mUueuee  ol  mjiiu- 
lieant  ami  unialli  iuii|naiitiliahle  miMiii;  pioeem-i 
I tun  it  ii  ilillieult  to  aiseis  the  i elat n e lole  ol 
ehenuitii  ami  ill  iiaunei  in  iletei  in  mine  the  ohn-t  t a 
non  ilata  \ln',  heeauie  ot  the  lonjiei  tune  leale. 
moie  mteileieuee  thiouith  eonieinoii  to  melt 
tempot  an  nnki  ,n  CIONO  . N O.  ami  1 1 NO  i shoulil 
he  eipeeteil  I lm  uieiean-il  mteileieuee  makei  it 
eien  mote  ilillieult  to  ioit  out  the  elieeti  ot  the 


Mli  Ift.  I » liu Ulnl  rt-i|'imir  a(  OII  mill  III'  In  * inlar  ei'liinr 
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catalytic  cycle  teacitonx  I malls,  at  sunrise  ami 
winxei  there  is  considerable  uncei taints  on  t Iso 
proper  treatment  lor  the  tlneet  and  Mattered  solar 
lluv  a>.  a tunction  ol  altitude  sshen  the  sun  i\  neat  the 
horizon  I hi»  then  leads  to  iirtccrtarritv  in  the 
representation  ot  photodissocuilion  processes  tn  the 
theotettcal  calculations  \t  least  tot  the  coming 
eclipse  this  will  not  he  a limitation  on  the  model, 
i e . the  eclipse  takes  place  neat  local  noontime 
We  would  like  to  teemphasi/e  the  point  that  the 
theotettcal  predictions  should  he  stewed  most  Is  in  a 
qualitative  sense  I he  rate's  ot  deer  ease  anil  increase 
in  mdistdual  ttace  species  concentrations  and  the 
relatise  changes  in  magnitudes  at  totalits  arc  the 
interesting  parameters  tot  comparison  in  the  meas 
tirements  It  an  espeiimental  piogtam  can  he 


estahlished.  and  sufficient  local  conditions  are 
knossn.  ness  calculations  can  he  carried  out  to  mils 
test  out  understanding  ot  stratospheric  chemistis 

Mthough  we  base  tocused  tnamls  on  the  local 
changes  at  mdistdual  altitudes,  death  total  column 
changes  are  also  observable  Since  the  variations  lot 
iiulisidual  species  at  all  altitudes  all  tollosv  the  same 
genetal  individual  trends,  the  total  column  lot  them 
also  s anes  m the  same  manner  I ot  example,  the  l It) 
column  shosvs  the  same  approximately  s nun  delas 
m reaching  minimum  as  is  shown  m I ig  IS.  and  the 
NO-  column  shosvs  a lactor  2 inctease at  totalits  l or 
those  interested  experimentalists,  a more  detailed 
report  including  x arious  changes  in  total  column  and 
mdixutual  sertical  ptofiles  is  asailahle 


3.  ASSKSSMKNT  STl  I)IKS 


3.1  Change  in  Ofttnr  because  of  SST 
Kmissions.  ( hlorofluoromethunes, 
and  Nitrous  Oxide 

I he  one-dimensional  transport  kinetics  model 
has  heen  used  to  assess  the  potential  chemical  cttects 
ol  sesetal  man-made  pettuihanons  to  the  stiatox- 
pherc  Whenexet  significant  changes  occur  to  the 
model  chemistis  or  the  treatment  ot  other  physical 
processes,  the  assessment  studies  are  repeated  in 
order  to  assess  the  eflect  ot  the  changes  on  the  model 
sensuous  1 he  results  presented  here  were  obtained 
using  the  latest  version  ot  the  model 

SSI  Kmivsions 

1 he  perturbation  xxe  use  as  a test  case  is  an  Nr), 
miection  at  a rate  ot  1000  molecules  cm  ' oxer  a I- 
km-thick  laser  centered  at  either  I'  ot  '0  km 
altitude  I his  is  equivalent  loan  annual  injection  rate 
ot  f>  2 \ 10'  kg  st  ot  NO  into  the  Northern  Hemis- 
phere l sing  the  NO,  emission  indexes  front  the 
Cl  \P  Report  ot  Findings  (tirobecker  <-r  <il . |0'4). 
this  iniectton  rate  is  approxtmatcls  equivalent  to  a 
tleet  ol  l2lH>  Concorde-txpe  SSI's  with  a enuse 
altitude  ot  1“  km  or  approximately  r>  H'  advanced 
SSI"'  with  a cruise  altitude  ol  '0  krn 

With  cuircnt  technologs  the  NO,  emission 
index  is  IS  g kg  ot  fuel  In  the  future  it  mas  be  pox 
si  Me  to  reduce  the  NO,  emission  index  to  h g kg  tuel 
In  addition  to  Nr),,  water  xapot  is  also  emitted  as  a 
product  ot  combustion  Sut  with  an  emission  index  ot 
I '50  g kg  tuel  Future  reductions  m the  Nr), 
emission  index  would  not  ctmnge  the  watei  xapot 


emission  index  Consequently,  tot  an  Nr),  miection 
rate  ot  1.1)00  molecules  cm  's  . the  H r)  miection 
tale  would  be  l".000  molecules  cm  x . tor  current 
teclmologs  and  531.000  molecules  for  tuture  tech- 
nology I he  latter  case  is  equivalent  to  increasing  the 
SS 1 tleet  xi/e  bs  a factor  ot  three  since  the  total  NO- 
miection  remains  unchanged 

1 he  computed  change  m total  o/one  is  given  in 
I able  3 lot  various  NO,  and  11-0  injections  into  an 
ambient  atmosphere  that  includes  cut  rent  levels  ot 
ehlorolluoromethanex  1 he  ambient  concentration 
ol  CIO,  is  I 10  ppbs 

Iniectmg  only  Nr),  results  in  an  increase  m total 
o/one  ot  slightls  more  than  21',  lot  both  I '-  and  20- 
km  uncct ions  these  increases  in  total  o/one  are 
larger  than  those  computed  a seat  ago  (I  tit  her  *•/,//.. 
Id?').  Iniectmg  ottls  water  xapot  into  the  stratos- 
phere results  in  a net  reduction  ot  total  o/one  In  the 

Table  2 Change  in  total  o/one  due  to  emissions  of 
NOv  and  HjO  at  injection  altitudes  of  I'  and  20  km 


Inio. non  rale 


Change  in  total  o/one. 


Specie* 

inoleenh-s 

em_-'»"* 

1 ’ km  injection 

20  km  injection 

Nt\ 

t .000 

l.dtl 

2 2 s 

HjO 

0 

Nl\ 

1 .000 

1.73 

1.82 

It  20 

1 77.000 

N**\ 

t .000 

t .VS 

0s*0 

HjO 

S.)  1.000 

calculations  in  which  both  NO,  and  H O aie  injected, 
the  etlcct  ol  the  water  \apoi  is  small compatcd  with 
that  of  NO,  I he  change  in  the  local  o/one  concen- 
tration is  shown  in  fig  I7  lor  the  I ’-kin  miection 
altitude 
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( 'hlorofluorome  thanes 

1 he  perturbation  we  considered  was  the  histor- 
ical release  ol  1’FM‘s  up  to  |d75  followed  b\  a 
constant  release  rate  at  ld'5  lesels  until  steads  state 
I he  steady-state  release  rates  are  assumed  to  be  H>N 
kt  vr  for  Cl  C'b  and  3S0  kt  \rlorCI  l'l  1 he  com- 
puted change  m total  o/one  at  steads  state  is  - I s d , , 
which  is  slightlx  greater  ,b  >n  the  o/one  reduction  we 
computed  a seat  agoll  uther  i t ul . Id”,  p 701  I he 
model  calculations  indicate  a present-das  reduction 
in  total  o/one  ol  7 7‘ , relatise  to  IdSp  1 he  change  in 
local  o/one  concentration  sersus  altitude  at  steads 
state  is  shown  in  Kig  IS 

Nitrous  Osidc 

Penitntication  resulting  trout  the  use  ol 
mitogen  base  leitili/ers  has  been  purposed  as  a 
source  ot  increased  N.l)  concentrations  in  the  atmos- 
phere (Crut/en.  |d’41  Increased  worldwide  use  ol 
nitrogen-base  fertilizers  mas  result  in  a doubling  ot 
VO  concentration  in  the  next  Ml- 100  sears,  but  the 
rate  of  met  ease  of  VO  is  highix  uncertain  We  base 
computed  the  change  in  total  o/one  lot  a range  of 
increases  of  the  VO  eoncentiation  in  the  tropo- 
sphere, and  the  icsults  are 'hossn  in  I ig  Id  \ small 
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Multiple  of  present  N ,0 

Ht..  I«  t h*nj[r  in  lout  o/one  itue  to  an  increase  in  VO  con- 
crnlralion  al  ihe  rarih's  surface 

increase  (v.0  V,  1 in  total  o/one  is  computed  for 
increases  in  N-O  up  to  a doubling  of  the  present 
tropospheric  concentration  (770  ppbsl  Further 
increases  in  N O result  in  a decrease  in  total  o/one 
(AO,  -7  0>  i tor  .1  \ N (1  and  At),  -7  41',  for  4 \ 
N 111 


3.2  Potential  Kffects  of  Space-Shuttle 
F missions 

In  Id'p,  sse  studied  the  potential  effeet  ot  IK'l  in 
the  exhaust  ol  the  space  shuttle  on  the  stratosphere 


Assuming  60  space-shuttle  flights  per  year,  we  cal- 
culated a 0.25%  reduction  in  total  ozone  at  steady 
state  using  the  LLL.  one-dimensional  transport- 
kinetics  model.  Subsequent  to  these  calculations, 
numerous  changes  occurred  in  the  measured  values 
of  chemical  reaction  rates  important  to  stratospheric 
chemistry.  Probably  the  most  important  of  these 
changes  was  the  measurement  of  the  rate  of  HO:  + 
NO  — OH  + NO;  by  Howard  and  Evenson  ( 1977). 
This  reaction  was  found  to  have  a rate  approximately 
20  times  faster  than  previously  estimated,  with  a 
measured  room  temperature  value  of  8 X 10  12 
cm’/sec.  In  addition  there  were  a number  of  smaller 
but  nonetheless  significant  changes  to  other  reaction 
rates  (due  to  new  measurements  or  estimates), 
particularly  to  reactions  involving  chlorine  species. 
Details  of  these  changes  can  be  found  in  Chang  et  at. 
(1977).  Because  of  these  changes  to  model  chemistry, 
it  is  useful  to  reexamine  the  potential  stratospheric 
effects  from  future  utilization  of  the  space  shuttle. 

Besides  emissions  of  approximately  60,000  kg  of 
HCI  per  flight  into  the  stratosphtre( based  on  NASA 
estimates),  there  are  also  other  exhaust  products  to 
be  considered.  There  are  expected  to  be  approxi- 
mately 260,000  kg  of  H;0,  130,000  kg  of  CO,  and 
4,200  kg  of  NO  emitted  per  flight  in  the  exhaust. 
Aluminum  oxides  emitted  in  the  exhaust  may  also  be 
impot  tant  to  particulate  formation,  but  the  impact  of 
this  effect  will  not  be  assessed  here.  Water  vapor  and 
carbon  monoxide  are  important  minor  constituents 
in  the  stratosrhere  with  mixing  ratios  by  volume  of 
approximately  4-5  ppm  for  H;0  and  30  ppb  for  CO. 
This  is  to  be  compared  with  approximately  1 .6  ppb  of 
total  odd  chlorine.  CIO,  (Cl,  CIO,  CIONO;,  and 
HCI),  calculated  for  the  present  atmosphere  in  the 
model.  Assuming  60  flights  per  year  for  the  shuttle, 
emissions  of  HCI  lead  to  approximately  a 5% 
increase  in  CIO*.  Nitrogen  oxide  emissions  are 
approximately  15  times  smaller  than  the  HCI  emis- 
sion, while  the  natural  NO,  (NO,  NO;  and  HNOt) 
concentrations  are  approximately  8-10  times  larger 
than  the  CIO,  concentration.  Therefore,  the  NO 
emissions  from  space  shuttles  should  be  insignificant 
relative  to  the  natural  concentrations.  The  emissions 
of  H;0  and  CO  from  space  shuttles  are  much  larger 
than  that  of  HCI,  but  their  ambient  concentrations 
are  also  much  larger  such  that  space  shuttles  should 
not  be  a significant  perturbation  to  water  vapor  and 
carbon  monoxide. 

Model  Resu\*s 

The  effect  or  chlorine  emissions  from  space 
shuttles  has  been  ‘ested  with  the  LLL  one- 
dimensional transport-kv,e(jcs  model.  The  chem- 
istry in  the  model  is  essential,;  ^at  recommended  by 
the  recent  evaluation  in  the  N AS*  rep0rt  on  Chloro- 


fluoromethunes  and  the  Stratosphere  (Hudson, 
1977).  The  hemispherically  averaged  HCI  source 
function  used  in  the  calculation  is  shown  in  Fig.  20. 
This  source  function  is  based  on  data  provided  by 
NASA  assuming  60  space  shuttle  flights  per  year. 
There  are  also  emissions  below  10  km,  but  these  are 
in  the  troposphere  where  it  is  assumed  that  the  HCI  is 
rapidly  removed  from  the  atmosphere  by  heterogen- 
eous or  rainout  processes. 


FIG.  20.  Hemispherically  averaged  HCI  emission  rale  in  the 
stratosphere  for  an  assumed  60  space-shuttle  flights  per  year. 


The  calculated  concentrations  of  atmospheric 
species  derived  for  current  levels  of  chlorofluoro- 
methanes(CFCL  and  CF;C1;)  were  used  as  the  initial 
conditions  for  this  study.  The  model  was  then  run  to 
steady  state  (approximately  20  years  of  model  time) 
using  the  assumed  space  shuttle  H Cl  emission  rate  in 
Fig.  20.  At  steady  state,  total  ozone  was  reduced 
0.27%  resulting  from  an  increase  in  atmospheric 
CIO,  of  0.08  ppbv  (a  5%  increase  in  stratospheric  odd 
chlorine).  The  calculated  change  in  local  ozone 
concentration  is  shown  in  Fig.  21.  The  maximum 
ozone  reduction  of  approximately  1%  occurs  near  40 
km.  This  is  probably  too  small  to  have  a significant 
impact  on  stratospheric  circulation  or  climate. 

The  calculated  change  in  total  ozone  is  very 
similar  to  that  computed  previously  in  1976.  Based 
on  calculations  of  the  effect  of  chlorofluoromethanes 
on  the  stratosphere,  a doubling  in  the  computed 
effect  of  ozone  would  have  been  expected  from  the 
change  in  the  rate  for  HO;  + NO  (Hudson,  1977). 
However,  other  changes  to  model  chemistry  since 
the  previous  calculations  have  compensated  for 
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from  \pm  o shuttle*  i'«Unlafion\ 

tins  In  particular,  reductions  in  the  rates  tor  the 
reactions  fit)  ) O * t 'I  * O- and  OH  t MCI  • H-O 
t C'l  have  occurred  1 he  rate  lor  CIO  t NO  • Cl  1 
NO-  has  increased  slightly,  and  the  N atom  sink 
(Ouewei  <7  ill . 1477)  tot  nitrogen  ovules  has  been 
1 educed  in  importance  due  to  rate  ehanges  In 
addition,  the  previous  calculations  had  a sniallei 
background  CIO,  concentration  because  present 
levels  ot  CI-Cli  and  Cl  -Cl  were  not  included  I he 
additional  ambient  chlorine  should  reduce  slight l\ 
the  computed  effects  ol  an  IICI  mieetion. 


J.J  Potential  Kffccts  of  Solar  Power 
Satellite  l aunch  Vehicles 

the  feasibility  ol  using  giant  solar  powei 
satellites  to  transmit  energy  to  the  earth  is  being 
studied  b\  the  Department  ol  Inergv  I he  system 
eurrcntlv  proposed  would  involve  as  many  as  1(H) 
solar  power  satellites  in  geosynchronous  mbit  l-'ach 
satellite  would  provide  5 ti\V  ol  electrical  powei 
which  would  be  beamed  to  receiving  stations  at  the 
ground  b\  microwave  transmission  Partial  assembly 
of  the  satellites  would  be  done  at  low-earth  orbit 
before  moving  them  to  a geosynchronous  orbit  A 


Heavy  i it t launch  Vehicle  (HI  I V)  would  be  used 
to  get  raw  materials  into  low-earth  orbit  Die  two- 
staged  (ot  biter  and  booster)  winged  launch  vehicle 
would  carry  a payload  of  about  4(H), 000  kg  Ihe 
booster  would  be  capable  ol  Hying  back  to  the 
launch  site  using  an -breathing  engines.  I he  01  biter 
would  glide  back  to  the  launch  site  upon  returning 
from  orbit.  CII,  and  (),-  would  be  used  as  fuel  for  the 
boostei  and  H - and  O-  would  be  used  lor  the  orbitet 
1 lie  pimiary  exhaust  products  ol  the  111  I V arc 
H -O  and  CO;  Some  NO,  would  be  produced  in  the 
trailing  shock  wave 

Using  preliminary  estimates  ol  the  injection 
rates  ol  NO,  and  H-O.  we  have  examined  the 
potential  effect  on  the  stratosphere  ol  HI  I V rocket- 
engine  emissions  I he  one-dimensional  transport- 
kinetics  model  was  used  lot  these  calculations  Using 
data  liom  Hoeing  lot  its  proposed  launch  vehicle  and 
assuming  141  llights  pet  year,  we  used  a source  rate 
lor  H-O  in  the  model  ot  I 4 \ |()4  metric  tons  km 
per  yeat  between  the  altitudes  ot  lb  and  SO  km  and 
b S S It)1  metric  tons  km  per  yeat  above  SO  km  Ml 
material  is  assumed  to  In-  deposited  into  the 
Northern  Hemisphere  in  the  calculation.  I he  source 
rate  for  CO-  should  be  comparable  to  that  lot  watei 
vapor,  but  this  is  expected  to  have  a negligible  effect 
on  the  amount  ot  t'O-  in  the  stratosphere  because 
ot  the  high  concentration  ot  carboiv  dioxide  leUvtisc 
to  that  lot  water  vapor.  Ihe  results  indicate  a 
maximum  increase  in  water  vapoi  ol  I1';,  whereas 
there  is  a lactoi  of  10  more  ambient  C’0>  than  H -O 
in  the  stratosphere 

I he  NO,  injection  rates,  shown  111  I ig  '2,  are 
based  upon  data  fiotn  N ASA  I angles  I stimatesare 
lot  solid  fuel  motors  ol  the  type  used  with  the  space 
shuttle  and  have  been  scaled  up  bv  the  nunibei  ot 
rocket  engines  to  be  used  in  the  III  I \ (lb)  and  bv 
the  i'lojectcd  nunibei  ot  llights  pci  yeat  (141)  I he 
NO,  inicction  rates  were  input  into  the  model  as 
nitric  oxide  Again,  the  injections  aie  assumed  to  In- 
in the  Northern  Hemisphere. 

Ihe  ambient  atmosphere  assumed  in  these 
calculations  is  that  containing  natural  sources  ol 
odd  chlorine  but  no  chloiotluoromcthanrs  I he 
added  chlorine  from  chloiolluoiomethanes  should 
have  little  impact  on  these  results  1 he  unpciturhcd 
concentrations  ol  H -O  and  NO,  are  given  in  I able  * 
Ihe  model  was  run  to  cuuilibrmm  with  the  source 
rates  described  above 

A change  m total  o/one  ol  -0  Old  "as  calcu- 
lated. I oca  I percent  changes  in  H-1*.  HO,,  and  Oi 
are  given  m I able  4 Ihe  n»«xiniuni  calculated 
change  in  local  o/one  was  < neat  the  top  ol  the 
stratosphere,  correspon*")K  •'*'  mciease  in  watei 
vapor  ol  I O V’i  I h."*ecrease  in  o/one  is  due  to  the 
added  HO,  resul»nl*  from  O('O)  t H O • 2011 
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N0X  injection  rate  — molecules/cm3  • s 


FIG.  22.  The  NO,  injection  rate  vs  altitude,  as  caused  by  39| 
HI.I.V  flights  per  year. 


Table  4.  Changes  in  species  concentrations. 


Altitude,  km 

Izxil  percent 

changes9 

ll2() 

HOxb 

NOxc 

03 

20 

0.20 

0.06 

0.04 

0 

25 

0.50 

0.46 

0.09 

-0.01 

30 

0.6: 

0.53 

0.07 

-0.01 

35 

0.74 

0.49 

0 

-0.10 

40 

0.84 

0.49 

0 

-O.I7 

45 

0.95 

0.48 

0 

-0.23 

50 

1.02 

0.44 

0 

-0.25 

55 

1.03 

0.41 

0 

-0.25 

aMany  of  these  percent  changes  are  based  on  the  last  signifi- 
cant figure  in  computer  printout  from  the  model. 

bHOx  = OH  + H02  + 2 X H202. 
cNOx  = NO  + N02  + HNO3. 


Perhaps  a more  important  effect  on  the  atmos- 
phere may  come  from  the  launch  vehicle  injections 
of  water  vapor  into  the  mesosphere,  since  this  might 
affect  ion  recombination  rates.  Because  our  model 
has  a maximum  altitude  of  55  km  and  docs  not  treat 
ion  chemistry  in  detail,  we  were  unable  to  examine 
this  effect. 


Table  3.  Unperturbed  species  profiles. 


Altitude, 

km 

Water  vapor 
mixing  ratio, 
ppmm 

NOx 

mixing  ratio, 
ppmv 

16 

2.28 

1.37  X I0~3 

18 

2.35 

2.59  X 10~3 

20 

2.45 

4.57  X 10~3 

22 

2.58 

7.40  X 10“ 3 

24 

2.71 

9.80  X I0-3 

26 

2.80 

1.13  X 10-2 

28 

2.88 

1.24  X I0“2 

30 

2.95 

1.34  X 10-2 

32 

3.01 

1.44  X I0-2 

35 

3.10 

1.57  X I0-2 

40 

3.23 

1.77  X 10~2 

45 

3.31 

1.82  X I0'2 

50 

3.35 

1.72  X 10-2 

55 

3.36 

1.64  x 10' 2 

3.4  Effects  of  a Massive  Pulse  Injection  of 
NO.  into  the  Stratosphere 

In  the  early  nineteen-seventies,  the  important 
role  of  NO,  in  regulating  stratospheric  ozone  was 
recognized  (Crutzen,  1970;  Johnston,  1971)  and 
potential  anthropogenic  sources  of  NO,  comparable 
to  the  natural  sources  were  identified.  These  anthro- 
pogenic sources  of  NO,  have  included  emissions 
from  supersonic  transports  (Grobecker  <7  al.,  1974; 
Johnston.  1971;  National  Research  Council,  1975a), 
N ;0  produced  as  a result  of  denitrification  from 
fertilized  soils  (Crutzen,  1974),  and  NO,  from  past 
nuclear  tests  (Foley  and  Ruderman.  1973)  or  from  a 
future  nuclear  war  (Hampson.  1974). 

Largely  as  a result  of  the  incorporation  of  recent 
measurements  of  chemical  reaction  rates  that  had 
previously  been  poorly  known  (Burrows  et  al.,  1977; 
DeMore  et  al.,  1977;  Howard  and  Evenson,  1978; 
Molina  et  al.,  1977).  model-predicted  ozone  reduc- 
tions from  NO,  sources  designed  to  simulate  a strat- 
ospheric aviation  effluent  (Broderick,  1977;  Popoff 
et  al..  1978;  Wucbblcs  et  al..  1977).  N;0  increases 
(Logan  et  al.,  1978)  or  nuclear  tests  (Chang  et  al.. 
1978)  have  been  greatly  reduced,  and  ozone  increases 
arc  calculated  for  some  scenarios.  However,  it  has 
been  estimated  that  the  NO,  input  to  the  strato- 
sphere after  a nuclear  exchange  might  be  nearly  two 
orders  of  magnitude  greater  than  that  produced  by 
the  other  hypothesized  NO,  sources  (MacCracken 
and  Chang,  1975;  National  Research  Council  1975a, 
b).  Because  the  processes  that  have  led  to  major 
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changes  in  model  sensitivity  to  other  types  of  NO, 
injection  might  he  overwhelmed  bv  such  perturba- 
tions, we  have  reinvestigated  model  response  to  very 
large  NO,  pulse  injections  comparable  to  those  that 
might  be  produced  in  a nuclear  exchange. 

We  have  used  our  one-dimensional  transport- 
kinetics  model  to  examine  i's  response  to  large 
pulse  NO,  injections.  I he  reaction  rates  are  the 
same  as  those  given  in  I able  A-l  (Appendix  A)  with 
the  exception  ol  the  reactions  listed  in  I able  5, 

1 he  perturbations  we  considered  consisted  ot 
pulse  injections. ol  NO,  into  the  layers  4 5-16  5 km 
(approximately  the  stabilization  altitude  range 
expected  lor  a 0.25-Mi  device).  1-  5-21  5 km  (the 
approximate  stabilization  altitude  range  expected 
for  a l-Mt  device),  and  17.5-28.5  km  (the  approxi- 
mate stabilization  altitude  range  expected  lor  a 
4-Mt  device).  It  was  assumed  that  6.7  N Id" 
molecules  of  NO  were  produced  per  megaton  ot  y ield 
(COM  IS  A.  1475)  and  that  the  NO  produced  is 
distributed  uniformly  over  the  Northern  Hemis- 
phere NO,  injections  corresponding  to  total  nuclear 
welds  ol  IDO.  l(H)t).  5000.  and  10000  megatons  were 
considered  for  each  of  the  three  injection  altitudes. 
1 hese  injection  scenarios  were  designed  to  explore 
the  effects  ol  a range  of  device  yields  and  total 
nuclear  yield  over  a range  that  might  be  plausible 
and  were  not  intended  as  realistic  scenarios. 


Results  and  Discussion 

I wo  computed  changes  m total  ozone  are 
shown  m I ig  ’.)  as  (unctions  ol  lime.  Very  large 
ozone  reductions  are  computed  loi  a total  y ield  ol  I (I4 
Ml  loi  I- and  4 Ml  devices  I hese  ozone  reductions 
dittei  little  Iront  those  we  obtained  m|475  loi  these 
device  yields  and  total  yield  ( Macl’i.icken  and 
C "hang.  1475)  In  conn  ast.  the  ozone  i eductions  com- 
puted toi  smallct  total  yield  oi  smaller  device  yield 
are  cuircntlv  much  smallci  than  computed  pre- 
V lotlslv 

I igmc  74  shows  the  computed  ozone  profiles  at 
0 05.  0 5.  and  I yeai  altei  the  inirctnm  loi  a 10*  Ml 
total  yield  ami  ditlerenl  device  yields  I lie  cones 
ponding  changes  in  local  ozone  concentration  are 
shown  in  I ig  7s  l ot  4 Mt  devices,  ozone  is  greatlv 
depleted  at  early  times  in  the  region  ol  the  initial 
miection.  and  substantial  changes  occm  in  the 
tioposphcie  as  a icsult  ol  the  increase  in  troposphei ic 
photodissociution  tales  even  belore  much  iniected 
NO,  is  tiansporled  into  the  troposphere.  At  later 
times  the  iniected  NO,  mixes  to  higher  and  lower 
altitudes,  reducing  ozone  at  high  altitudes  while 
increasing  ozone  m the  lowei  stratosphere  and 
troposphere  Sunilai  cl  led  s are  show  n for  the  results 
using  l-Mt  devices,  although  tropospheric  ozone  is 
allccted  bv  NO,  even  at  0.05  years  1 or  0.25-Mt 
devises  much  of  the  perturbation  appears  m the 


l uble  5.  Reaction  rates  that  differ  from  1**77  Chemistry . 
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FIG.  23.  Calculated  change  in  total  ozone  vs  time  for  various 
total  nuclear  yields  in  megatons.  Calculations  apply  to  individual 
device  yields  of  0.25,  1.0.  and  4.0  megatons. 

troposphere  immediately  after  the  injection,  and 
relatively  small  changes  in  o/one  occur  above  the 
height  of  injection. 

To  understand  the  small  change  since  1975  in 
model  sensitivity  to  very  large  NO,  injections  in  the 
middle  stratosphere,  it  is  necessary  to  understand  the 
sources  of  altered  model  sensitivity  to  smaller  injec- 
tions. Two  factors  have  been  materially  affected  by 


FIG.  24.  Calculated  o/one  concentration  vs  altitude  after  a 10* 
Mt  total  yield  pulse  for  various  device  yields. 

changes  in  chemical  reaction  rates  involving  HO; 
used  in  the  model. 

Use  of  3 X 10  " cm'  sec  (Dc More  el  al.,  1977) 
rather  than  2 X 10  cm'  sec  for  the  reaction 

HO:  + OH  - H O + O 

has  substantially  increased  the  computed  concentra- 
tions of  HO,  radicals  in  the  model  stratosphere.  Asa 
result,  reactions  of  HO,  radicals  with  O and  Oi  are 
more  important  sinks  for  o/one  in  current  models 
than  they  used  to  be. 

Use  of  8 X 10  1 cm'  sec  (Howard  and  Evenson, 
1978)  rather  than  2 X 10  " cm'  sec  for  the  reaction 

HO:  + NO  - NO:  + OH 
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has  caused  lit),  and  NO,  eataly /ed-o/onc  destine 
tn'ii  cycles  m interfere  with  cadi  othei  to  a grcatci 
extent,  ami  has  increased  the  ellicienci  ol  i lie  \(>, 
enhanced  nicthane"smog"  nicchanistn  lot  pioducing 
"/one  both  HO,  catuly/cd-o/onc  destruction  and 
methane  "smog"  o/one  production  are  important  in 
the  lowet  sitalospherc,  and  the  use  ol  nenei  data  lot 
lilt  reactions  has  led  to  a situation  in  which  ealeu 
la  ted  o/one  increases  heU'w  about  70  Km  can  be 
conipaiahle  to  o/one  decteases  computed  in  the  70 
45  Km  legion  (NO,  induced  o/one  changes  abuse 
45  Km  ate  small  loi  NO,  uiiectious  below  >0  Km. I 
I lieie  is  also  an  inteilerence  between  NO.  and  CIO, 
o/one  destruction  cycles  as  a result  ol  the  reactions 

NO  t CIO  • NO  t Cl 


NO  * CIO  I M • CINOi  t M 

It  should  be  noted  that  o/one  pi  oduction  \ i.i  the 
"smog"  mechanism  is  ultimately  limited  by  the 
methane  Iluv.  and  it  o/one  destruction  s in  HO,  oi 
l IO,  is  siippiessed  to  the  point  u becomes  negligible. 
Im  t het  suppression  ot  t liesc  processes  w ill  not  cause 
tut l het  D>  increases  thus,  the  NO,  induced 
increases  in  lowet  stratospheric  o/one  cannot 
be  increased  without  limit  Out  largest  injection 
scenarios  result  in  lOO-fold  increases  in  the  stratos- 
pheric burden  ol  NO.  and  more  than  I (Mold 
incieases  m the  burden  ol  NO,  ( NO  t NO  t UNO,  t 
NO,  t 7N  O.V  It  a large  NO,  pulse  reaches  the 
middle  stratosphere,  it  will  cause  a very  substantial 
computed  o/one  reduction  which  can  overwhelm  the 


relatively  modest  o/one  increases  computed  in  the 
lowet  stiatospheic  and  tioposphere 

ll  NO,  is  miected  at  a low  enough  altitude,  most 
ol  it  is  t ranspoi led  into  the  tioposphere  pat tlv  in  the 
lot m ol  UNO.  and  NO  whete  it  is  assumed  to  be 
removed  In  raiiiout  ptocesses  Vs  a result,  veiv  little 
is  Iranspottcd  into  the  midstintospherc  where  it 
would  ellectivelv  reduce  o/one  I Inis  the  computed 
o/one  incieases  in  the  lowet  stratosphete  and  tropo- 
sphere are  the  dominant  elicits  even  lot  vetv  large 
NO,  injections  below  appioximntclv  I 5 Km 

( oneliisioiis 

I i win  these  results  we  see  that  in  contiast  to 
othet  perturbation  scenarios  vv  e continue  to  compute 
veiv  large  reductions  in  stiatosplieric  o/one  tollovv 
mg  a massive  pulse  injection  ot  NO,  aKin  to  that 
expected  it  several  thousand  megatons  ot  nucleai 
explosives  m the  megaton  range  were  to  be 
detonated.  Ilovvcvet.  stnallei  uiiectious  ot  imeetions 
I tom  devices  ot  0.7.5  \u  oi  less  yield  would  be 
expected  to  have  a stnallei  etlecl  on  the  total  column 
ol  stratospheric  o/one 

I lie  o/one  column  recovers  more  rapidly  in 
these  calculations  than  it  did  in  om  eaihei  calctila 
tions  ( Macl'racKeu  and  ( hang.  Id’s)  | ,,  tl  Nt-|V 
large  extent  this  is  a retied  ton  ol  out  use  ol  stibstan 
ttally  tastei  tianspoit  above  77  Km  than  in  thecat  liei 
ealeulatioiis  l lie  giowth  ol  lovvei  stratospheric 
o/one  increases  is  a lesset  eontnbutoi  to  the  tastei 
recoverv 

V maioi  limitation  ,'t  these  calculations  is  that 
no  attention  is  paid  to  itUeiaciions  ol  changes  m 
stratospheric  chemistiv  with  the  atmospheric  simc 
lute  l ,'i  perturbations  as  l.uge  ,is  some  ,'t  those 


I 
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calculated  here,  such  interactions  might  be  impor- 
tant. Because  of  this  as  well  as  the  other  limitations  of 
one-dimensional  models  (Chang  el  til.,  1 978).  the 
computed  recovery  rates  should  he  interpreted  cau- 
tiously. I'he  computed  peak  ozone  depletions  are 
only  weakly  dependent  on  the  transport  parameters. 
K,.  for  the  very  large  injection  scenarios  at  the  higher 
altitudes.  However,  both  the  recovery  time  and  the 
computed  sensitivity  to  intermediate  injections  are 
sensitive  to  K.,  since  K,  influences  the  residence  time 
of  in  jected  NO,  and  the  ambient  distributions  of  NO, 
and  other  species. 

3.5  Effects  of  Stratospheric 

Perturbations  on  the  Earth's 
Radiation  Budget,  and  the  Climatic 
Implications 

Since  the  end  of  the  CIAP  program  in  1975. 
there  have  been  many  changes  to  the  numerical 
models  used  to  compute  the  effects  of  stratospheric 
perturbations  on  0«.  New  species  (mostly  chlorine 
containing)  and  new  reactions  have  been  added. 
Poorly  known  reaction  rates  have  been  remeasured 
using  more  accurate  methods,  and  the  Oi  photo- 
dissociation rate  calculation  has  been  improved. 
Temperature  coupling  and  feedback  on  reaction 
rates  have  been  added,  and  diurnal  averaging  has 
been  accounted  for  more  accurately.  These  changes 
have  significantly  modified  the  Oi  impact  assess- 
ments. It  is  important,  therefore,  to  re-analv/e  the 
impact  of  the  potential  changes  in  stratospheric 
composition  on  the  global  radiation  budget. 

In  a previous  assessment  of  the  effects  of 
changes  in  Oi  and  NO;  concentrations  on  solar 
absorption  and  stratospheric  heating  rates,  l.uther 
(1976)  used  species  concentration  profiles  derived 
from  our  then-current  transport-kinetics  model.  In 
updating  the  assessment,  it  is  desirable  to  present  the 
results  in  a more  general  way  so  that  the  radiative 
effects  can  be  inferred  for  any  transport-kinetics 
model  results.  This  would  also  enable  the  radiative 
effects  to  be  upgraded  simply  whenever  there  is  a 
significant  change  in  the  stratospheric  model  assess- 
ments. We  have  attempted  to  do  this  for  calculations 
of  total  solar  absorption  by  Oi  and  NO;. 

Comparison  of  Changes  in  Solar  Absorption  by 
O.  and  NO. 

In  the  case  of  a stratosphere  perturbed  by 
an  NO,  injection,  l.uther  (1976)  showed  that  the 
increase  in  solar  absorption  by  NO;  at  steady  state 
was  a significant  fraction  (35  to  51K< ) of  the  decrease 
in  solar  absorption  by  Oi.  Since  that  time,  the  sen- 


sitivity of  Oi  to  an  NO,  injection  has  decreased  in  the 
transport-kinetics  models.  Consequently,  the  change 
in  solar  absorption  by  NO;  is  now  expected  to  be  a 
much  larger  fraction  of  the  change  in  solar  absorp- 
tion by  Oi. 

Solar  absorption  by  Oi  is  show  n in  Tig.  26  as  a 
function  of  0«  column  density.  The  absorption  rate 
given  is  the  instantaneous  \alue  for  a solar  zenith 


Total  ozone  — atm  • cm 


FIG.  26.  Solar  absorption  by  o/one  for  a solar  zenith  angle  of 
60  A,  is  the  surface  albedo. 
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FIG.  27.  Solar  absorption  by  NO  for  a solar  zenith  angle  of 
60". 


angle  ol  60  I he  radiative  transfer  model  used  to 
compute  the  solar  absorption  rate  includes  Rayleigh 
scattering  and  assumes  a cloudless,  plane-parallel 
atmosphere  above  an  isotropically  scattering  ground 
(luther.  1076)  I he  solar  absorption  by  NO  is 
shown  in  I ig  27  also  tor  a solar /enith  angle  of  60' 
Ambient  and  perturbed  species  concentration 
profiles  were  computed  using  the  1 I I one- 
dimensional transport-kinetics  model  I he  model 
chemistry  is  listed  in  I able  A-l.  except  for  a few 
reactions.  I he  rates  used  lot  the  reactions  Oi  + HO  , 
Oi  + H.  HO:  + HO  . and  Cl  * HO  - are  given  in  I able 
A-2.  I he  rate  ol  0.0  A It)  was  used  for  the  reaction 
NO  + HO.,  and  a rate  of  5.0  X 10  " exp(-l238  I > 
was  used  for  CIO  + CIO.  I he  reactions  ()(lD)  + Oi 
and  HO;  -1-  HO:  + H;0  were  not  included.  I wo  per- 
turbation eases  were  considered:  NO,  injections  at  17 
or  20  km  at  the  rate  of  1000  molecules  cm  s 1 uni- 
formly distributed  over  a l-km-thick  layer.  This 
corresponds  to  a hemispheric  injection  rate  of  6.2  X 
10"  kg  yr  of  NO:.  Changes  in  the  O.  and  NO:  column 
densities  and  solar  absorption  rates  for  a solar  /enith 
angle  of  60°  and  a surface  albedo  of  0.25  are  sum- 
marized in  Table  6.  The  unperturbed  column  densi- 
ties are  if. 730  X 10  ' molecules  cm  (0.382  atni.em) 
for  Oi  and  3.506  X It)1  molecules  cm  tor  NO:. 


Table  6.  Increases  in  total  atmospheric  solar  absorp- 
tion by  Oy  and  NOs  due  to  NOx  injections  at  the  rate 
of  6.2  X 108  kg/vr.  Calculations  are  for  a solar  zenith 
angle  of  60'  and  a surface  albedo  of  0.25. 


Quantity 

NOv  injection  altitude 

1 7 km 

20  km 

at)  3 

1.06% 

0.39% 

aNOj 

6.87% 

1 5 .60'i 

AAhsjO  j) 

0.16  Vi/m2 

0.06  W m‘ 

AAbs(N02) 

0.06  W/m2 

0.1  JW  m' 

In  both  eases  the  injection  ol  NO,  resulted  in  a 
small  increase  in  total  Oi.  Although  the  changes  in 
the  Oi  column  density  are  small,  the  change  in  the 
local  Oi  concentration  (Tig.  28)  may  be  several 
percent.  I he  change  in  column  density  relleets  the 
net  difference  between  regions  ol  Oi  increase  (below 
22  km)  and  Oi  decrease  (above  22  km)  Conse- 
quently. although  the  change  in  the  net  heating  may 
be  small,  the  redistribution  in  altitude  of  where  this 
heating  occurs  may  be  significant.  I he  change  in  the 
local  concentration  ol  NO  is  shown  in  I ig  2d 


Change  in  ozone  concentration  — % 

til. , ’s.  Change  in  local  o/one  concentration  due  to  NO. 
injections  (as  NO  ) of  \ III'  kg  ir. 


Hit.  I1*.  Change  in  local  NO  concentration  due  to  NO.  injec- 
tions (as  NO  t of  6.2  X 10'  kg  rr. 


Changes  to  the  Solar  Radiation  Budget  and  Heating 
Kates 

In  addition  to  the  two  NO,  injection  eases 
already  described,  we  also  considered  the  steady- 
state  Oi  reduction  due  to  Cl  M's  at  the  Id 75  release 
rate.  1 he  change  in  O:  column  density  computed  at 
steady  state  was  -14. 4‘  < l or  a solar  /enith  angle  ol 
60  and  a surface  albedo  of  0.25.  there  is  a reduction 


in  solar  absorption  b>  0«  of  2.45  W nr  1 he  change 
in  the  local  O.  concentration  is  shown  in  f ig  30. 

1 here  is  a reduction  m 0>  concentration  at  all  alti- 
t uiles  w ith  the  largest  percent  reduction  occut  ring  at 
40  km 

In  addition  to  affecting  the  gaseous  absorption 
ol  solar  radiation  in  the  stratosphere,  changes  in 
stratospheric  composition  also  affect  solar  absorp- 
tion m the  troposphere  by  changing  the  atmospheric 
transmissivity  A decrease  in  O.  column  density . for 
example,  allows  more  solar  radiation  to  reach  the 
troposphere,  thereby  increasing  both  the  amount  of 
solar  absorption  and  the  amount  of  radiation  scat- 
tered back  to  space.  Changes  to  the  solar  radiation 
budget  are  presented  in  I able  7 for  the  three  pertur- 
bation cases  I lie  net  solar  flux  at  the  top  of  the 
atmosphere  is  defined  as  the  incoming  flux  minus  the 
outgoing  flux.  I he  values  for  the  change  in  solar 
absorption  in  the  troposphere  include  the  change  in 
gaseous  absorption  and  the  change  in  absorption  by 
the  earth's  surface.  I he  changes  in  the  instantaneous 


Hi;.  HI  l hangr  in  oronr  concentration  «t  straits  stair  tine  In 
( KM\  at  the  l*>7S  release  raw 


solar  absorption  rates  for  the  NO,  injection  eases  are 
small  compared  with  the  unperturbed  values.  1 he 
t'l  M perturbation  ease,  on  the  other  hand,  repre- 
sents a significant  perturbation  to  the  solar  radiation 
budget.  I here  is  a significant  reduction  in  the  solar 
absorption  in  the  stratosphere  and  an  increase  in  the 
troposphere. 

I he  changes  in  the  Ot  and  NO-  solar  heating 
rates  arc  presented  in  l;ig.  31  for  the  17-km  NO, 
injection  and  in  f ig.  32  for  the  20-km  NO,  injection, 
t he  Oi  and  NO;  both  contribute  to  an  increase  in  the 
heating  rate  in  the  lower  stratosphere.  Consequently, 
the  change  in  solar  heating  tends  to  increase  the  tem- 
perature in  the  upper  stratosphere. 

1 he  change  in  the  solar  heating  rate  due  to  a 
CIM  perturbation  is  shown  m f ig.  33.  1 he  reduction 
m the  solar  healing  rate  in  the  lower  stratosphere  is 
approximately  the  same  magnitude  as  the  increase  m 
heating  rate  computed  for  the  NO,  injection  cases 
I he  change  in  the  solar  heating  rate  is  -0.07  k day  at 
24  km.  -0.02  k day  at  20  km.  and  -0.01  k day  at  IS 


Change  in  solar  heating  rate  — K day 


llli  M . i hange  in  (he  insdntiwtw  volar  tinting  ratev  (ot  * 
volar  zenith  angle  of  f»fl  at  \ 0,2'  due  to  an  NO,  injection 

at  P km 


fable  7.  The  perturbed  solar  radiation  budget  (solar  zenith  angle  r-  t>0  , surface  albedo  = 0 251 


Perturbation  case 

ft.  1 

Quantity 

Nt\.  1 7 km 

NOv  20  km 

CIM 

l nper  turned 
value 

.ANct  solar  flux  at  lop  of  atmosphere 

0.1  W m’ 

0.0  Vi  mJ 

-1.0  W nC 

470  W 

•XSttMosphctic  *>lar  absorption 

0.13 

o.os 

-2.35 

34.1 

•X  Troposphere  surface  .slat  ahmrption 

-0.1 

-O.l 

1.3 

43b 

Change  in  solar  heating  rate  — K/day 


H(.  '2.  ( han{r  in  (hr  instantaneous  solar  heating  rater  for  a 
solar  zenith  angle  of  (>0  and  \ (1.25  due  to  an  VO.  injection. 


H(i.  t hange  in  the  solar  healing  rate  of  ozone  at  steads  state 
for  a solar  zenith  angle  of  hi)  ,td  V 0.25  due  to  t't  M's  at  the 
IO15  release  rate. 

knt  In  this  ease  there  is  cooling  at  all  altitudes  with 
the  largest  effect  in  the  upper  stratosphere 

Potential  ( lunatic  Impart 

t he  climatic  impact  ol  changes  in  stratospheric 
composition  depends  upon  hs>th  the  solar  and  long- 
wave effects  ol  the  perturbation.  Changes  111 
temperature,  in  addition  to  affecting  the  transfer  of 
longwave  radiation,  also  affect  atmospheric  stability 
and  transport.  Here  we  will  consider  only  the  global 
impact  of  the  solar  and  longwave  effects  and  neglect 
am  potential  feedback  on  transport. 


A previous  assessment  of  the  effect  ol  changes  in 
Oi  and  NO-  on  surface  temperature  bv  Ramanathan 
<7  a!  ( I07M  showed  a cooling  at  the  surface  associ- 
ated with  a reduction  in  O.  due  to  an  NO,  iniection 
Only  reductions  in  O;  were  considered  (not  111 
creases),  and  it  was  assumed  that  the  changes  m O 
and  NO  were  uniform  (percentagewise)  between  17 
and  -JO  km.  Our  present  results  differ  tiom  these 
modeling  assumptions  in  that  the  changes  in  O.  and 
NO-  concentrations  are  not  uniform  with  altitude, 
and  there  is  a net  increase  in  Oi  column  density  lailtet 
than  a decrease.  Nevertheless,  the  work  of  Ramana- 
than  (7  id.  (|vJ7(v)  is  useful  because  it  demonstrates 
the  importance  ol  the  longwave  effect  ol  the  per- 
turbation. I he  reduction  in  total  O,  tended  to  warm 
the  troposphere  by  increasing  the  transmissiv itv  of 
the  stratosphere  for  solar  radiation.  I he  reduction  in 
stratospheric  temperature  due  to  reduced  Oi. 
however,  had  a greater  effect  on  the  longwave 
radiation  emitted  downward  from  the  stratosphere. 

I he  net  result  was  a slight  reduction  in  surface 
temperature. 

A similar  calculation  using  our  current  model- 
ing results  would  predict  an  increase  in  temperature 
in  the  lovvei  stratosphcrcfl  utherand  Duevver.  1 078 f. 
It  is  possible  that  the  longwave  effect  would  also 
dominate  in  this  case.  Although  the  net  change  in 
temperature  is  uncertain,  it  is  clear  that  the  longwave 
ellect  would  tend  to  warm  the  troposphere,  whereas 
the  solar  effect  would  tend  to  cool. 

In  attempting  to  assess  the  effect  of  changes  in 
stratospheric  composition  on  the  global  mean 
surface  temperature,  it  is  desirable  that  the  same 
model  be  used  throughout  the  study  Since  vv c do  not 
have  a climate  model  that  may  be  applied  to  this 
studv . we  will  rely  on  the  results  of  other  researchers 
I he  two  methods  used  here  are  the  published 
results  iif  Kamanathan's  radiative  equilibrium  model 
and  Budyko's  (Idl'd)  empirical  formulation  relating 
the  change  in  outgoing  longwave  flux  al  the  top  of 
the  atmosphere  to  the  change  in  surface  temperature. 

According  to  Kamanathan's  model  (Ramana- 
than <7  ill..  I '*"’(').  the  change  in  surface  temperature 
is  related  to  the  change  in  o/one  bv 

A I (tVOOd  to  0 014  kt  AO,  . 

where  AOi  is  the  pet  cent  change  in  o/one  concentra- 
tion applied  uniformly  over  the  altitude  region  1 2-40 
km  l ire  first  coefficient  applies  to  the  assumption  ol 
constant  cloud  top  altitude  whereas  the  second 
coefficient  applies  to  the  assumption  of  constant 
cloud  top  temperature.  A similar  expression  relates 
the  change  in  surface  temperature  to  the  percent 
change  in  NO.  concentration  between  17  and  40  km 


AT,  = (3  X 10  4 to  6 X 10  4 K)  ANO:  . 


Table  9.  Estimated  change  in  surface  temperature 
using  different  assessment  methods. 


The  expressions  are  not  applicable  directly  to  our 
results,  but  they  may  be  used  to  estimate  an  approxi- 
mate  upper  limit  on  AT,. 

The  present  changes  in  O,  and  NO:  averaged 
over  various  altitude  ranges  between  12  and  40  km 
are  given  in  Table  8.  Since  the  change  in  species  con- 
centrations in  the  lower  stratosphere  will  have  a 
greater  effect  on  the  downward  longwave  flux  into 
the  troposphere  than  changes  at  higher  altitudes,  we 
can  get  a rough  estimate  of  AT,  by  using  the  values  of 
AOi  and  ANO:  for  the  12-22  km  region.  Using  these 
values  tends  to  overestimate  AT,  because  it  assumes 
that  the  ozone  concentration  has  increased  by  the 
same  percentage  at  higher  altitudes  when  the  calcula- 
tion indicates  that  smaller  changes  occur  at  higher 
altitudes. 

Table  8.  Percent  change  in  O3  and  NO  2 concentra- 
tions averaged  over  various  altitude  ranges. 


Average  over 
altitude  range,  km 


Perturbation 

Quantity 

12-22 

22-40 

12-40 

NOx  injection  at  17  km 
at  rate  of  6.2  X 10® 

AO3,  % 

4.3 

-0.7 

0.7 

kg/yr 

AN02,% 

30 

6.0 

8.7 

NOx  injection  at  20  km 

AO3.  % 

4.9 

-2.0 

0.0 

at  rate  of  6.2  X 10® 
kg/yr 

aN02.% 

56 

15.6 

20.1 

Values  of  AT,  computed  using  the  relations 
above  and  values  of  AO>  and  A NO:  forthe  12-22  km 
region  are  given  in  Table  9.  The  change  in  surface 
temperature  is  estimated  to  be  less  than  0.1  K for 
both  the  1 7- and  20-km  injection  altitudes.  For  com- 
parison. values  arc  also  given  in  fable  4 for  AT,  using 
Budyko’s  empirical  formulation 

AT,  = 0.701  K AF  , 

where  AF  is  the  change  in  the  outgoing  longwave 
flux  in  W m at  the  top  of  the  atmosphere.  The  value 
for  AF  is  obtained  by  dividing  in  half  the  change  in 
net  solar  flux  at  the  top  of  the  atmosphere  given  in 
Table  7 to  account  for  day-night  averaging.  At  equil- 
ibrium the  change  in  the  outgoing  longwave  flux 
equals  the  change  in  the  net  incoming  solar  flux. 
Again  the  computed  values  of  AT,  are  significantly 
less  than  0. 1 K. 

The  same  methods  were  applied  to  the  CFM 
perturbation  results.  In  this  case,  however,  the  value 


Change  in  surface 
temperature,  K 

Perturbation 

Ramanathan 

Budyko3 

NOx  injection  at  1 7 km 

AO3 

ANO2 

AO3  + ANO2 

0.04  to  0.06b 
0.01  to  0.02b 

0.04 

NOx  injection  at  20  km 

AO3 

ANO2 

AO3  + ANO2 

0.04  to  0.07b 
0.02  to  0.03b 

0.02 

CFM’sat  1975  release  rate 
AO3 

ACFCI3.  ACF2Cl2C 

-0.13  to  -0.21 

0.44 

-0.35 

0.44 

aATs  = 0.701  K AF(W/m2). 

'’The  method  used  to  determine  these  values  tends  to  over- 
estimate the  change  in  surface  temperature. 


cFrom  Wang  el  al.  (1976).  ATS  = 0.38  to  0.57  K. 


of  AOi  averaged  over  12-40  km  (Table  8)  was  used  to 
estimate  AT,.  The  resulting  values  of  AT,  are  pre- 
sented in  Table  9.  I nt+iis  case  the  change  in  the  global 
mean  temperature  is  on  the  order  of  sev  eral  tenths  of 
a degree,  w hich  may  be  a significant  clima  te  pertur- 
bation since  regional  changes  in  temperature  might 
be  several  times  greater  than  the  global  mean. 
Pollack  el  al.  ( 1976)  suggest  that  0. 1 K is  a threshold 
value  for  considering  the  change  in  global  mean  sur- 
face temperature  to  be  significant.  Computing  v alues 
of  AT,  associated  with  past  major  changes  in  clim- 
ate. they  infer  that  values  below  0. 1 K would  not  have 
major  consequences. 

In  the  case  of  a CFM  perturbation,  there  is  the 
additional  effect  of  the  longwave  properties  of 
CF:CI'  and  CFCIi.  At  steady  state  we  compute  the 
tropospheric  concentrations  of  the  CFM's  to  be 
larger  than  present  day  levels  by  0.6  ppbv  for  CFCIs 
and  1 .52  ppbv  for  CF:C1;.  The  greenhouse  effect  of 
CFM’s  has  been  studied  by  Ramanathan  ( 1975).  and 
he  finds  the  following  relationships  between  the 
change  in  surface  temperature  and  the  change  in 
tropospheric  CFM  concentration: 

AT,  = 0.215  K ACF-CI;  (ppmv) 

AT,  = 0.180  K ACFCI.  (ppmv)  . 
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Ramanathan’s  model  results  therefore  imply  an 
increase  in  surface  temperature  of  0.44  k due  to 


! 


" ■ 


C'KM's.  which  is  the  same  magnitude  as  predicted 
using  Budyko's  empirical  formulation.  A similar 
calculation  has  been  performed  by  Wang  ctol.  (1976) 
in  which  an  increase  in  surface  temperature  ol  0.38  to 
0.57  k due  to  these  t'FM  concentrations  is  predicted. 

1 his  result  brackets  the  value  ol  0.44  k obtained  by 
the  other  methods. 

We  now  consider  the  potential  climatic  ettect  ol 
the  water  vapor  injected  into  the  stratosphere  along 
with  the  NO,  from  SS  I engine  emissions.  Assuming 
an  emission  index  of  1250  g kg  fuel  for  water  vapor 
and  18  g kg  fuel  for  NO,,  the  associated  water  vapor 
injection  rate  would  be  4.3  x [01  kg  yr.  According 
to  our  model  calculations  in  which  the  stratospheric 
water  vapor  profile  is  computed,  this  would  lead  to 
an  increase  in  the  stratospheric  water  vapor  mixing 
ratio  ot  0.1  ppmtti  for  a 17-kni  injection  and  an 
increase  of  0.31  ppmm  for  a 20-km  injection. 
Kamanathan's  model  was  also  used  in  the  CIAP 
study  to  estimate  the  change  in  surface  temperature 
resulting  from  a change  in  stratospheric  water  vapor 
mixing  ratio  (Grobecker  ci  ol..  1974): 

Al.  (0.2  to  0.3  k)  AH.O  (ppmm)  . 

I he  estimated  change  in  surface  temperature  due  to 
these  changes  in  stratospheric  water  vapor  are  given 
in  fable  10  fhese  values  of  AT,  are  larger  than  those 
due  to  AO.  and  ANO  resulting  from  SSI  engine 
emissions,  but  they  are  still  less  than  0.1  k. 

In  addition  to  NO,  and  water  vapor.  SS'I 
engines  also  emit  SO.,  which  is  converted  to  sulfate 
aerosols.  Assuming  an  emission  index  of  1.0  g kg 
fuel,  the  SO  injection  rate  would  be  3.4  X 10  kg  yr 
when  the  NO,  injection  rate  is  6.2  X 10*  kg  yr.  fhe 
change  in  surface  temperature  due  to  increased  stra- 
tospheric aerosols  ( 75‘ , II  SO»l  has  been  computed 
h\  Pollack  cl  ol.  (1976)  t sing  a radiative-convective 
model,  they  find 

AT.  = (-6.3  to  -It)  k)  Ar  . 


Table  10.  Change  in  surface  temperature  due  to  an 
increase  in  stratospheric  water  vapor. 


Altitude  of 

UjO  in 

injection. 

stratosphere. 

A f s* 

km 

ppmm 

K 

17 

0.1 1 

0.02  -0.03 

20 

0.31 

0.06-0.09 

where  Ar  is  the  increase  in  the  stratospheric  optical 
depth.  Pollack  ci  ol.  find  that  Ar  is  related  to  the 
mass  density  (m)  of  surfatc  aerosols  in  *ig  m'  aver- 
aged over  a 10-km  thick  layer  by  the  expression  Ar  = 
0.03 1 nr.  the  expression  Ar  = 0.038  m was  used  in  the 
CIAP  Report  of  Findings  (Grobecker  ft  al . 1074). 
The  aerosol  mass  density  is  obtained  from  the 
expression 

m — HFtc  M,  M,  . 

w here  H is  the  fraction  of  the  aerosols  deposited  in  a 
given  hemisphere.  F is  the  SO:  emission  rate,  t is  the 
residence  time,  c is  the  conversion  efficiency.  M.,  is 
the  molecular  weight  of  the  aerosol,  and  M,  is  the 
molecular  w eight  of  SO:.  For  a sulfuric  acid  solution 
that  is  75‘;  H.SOj  by  weight.  M.  M,  = 1.6.  Values 
for  the  various  quantities  used  by  Pollack  cl  al. 

( 1976)  and  used  in  the  CIAP  Report  of  Findings  arc- 
shown  in  fable  1 1,  fhe  resulting  values  for  Ar  are 
considerably  smaller  (by  nearly  a factor  of  3)  using 
Pollack  ci  al.'s  values.  In  both  cases  the  values  lor 
AT,  are  estimated  to  be  less  than  -0.01  k for  a 1 7-ktn 
injection  altitude  and  less  than  -0.02  for  a 20-km 
injection  altitude. 

Conclusions 

Changes  m O,  and  NO:  concentrations  due  to 
an  NO,  injection  lead  to  increases  in  the  total  solar 
absorption  of  both  species.  I hese  increases  are  small, 
however,  when  compared  to  the  total  energy 
absorbed  by  the  stratosphere,  and  the  climatic  effect 
is  estimated  to  be  negligible.  I sing  the  criterion  that 
changes  in  surface  temperature  less  thanO.  I k would 
not  have  major  consequences,  none  of  the  SSI 
engine  emissions  (NO,.  H.-O.  and  SO;)  are  estimated 
to  have  a major  climatic  etfect.  fhe  largest  individual 
effect  on  surface  temperature  is  that  ol  water  vapor, 
w hich  is  estimated  to  cause  a temperature  increase  of 
0.06-0.09  k for  a 20-km  injection  at  4.3  X 101"  kg  yr. 
The  combined  effect  of  all  engine  emissions  on 
climate  is  likely  to  be  an  increase  in  global  mean 
temperature  of  less  than  0.1  k. 

In  the  case  of  a CFM  perturbation,  there  is  a 
significant  impact  on  the  solar  radiation  budget 
leading  to  a decrease  in  the  solar  heating  of  the  stra- 
tosphere and  an  increase  in  solar  heating  of  the 
troposphere  and  earth's  surface.  I here  is  also  an 
increase  in  the  solar  flux  scattered  back  to  space, 
fhere  is  considerable  uncertainty  about  the  change 
in  surface  temperature  at  steady  state  caused  by  the 
reduction  of  Oi  due  to  CFM's  released  at  the  1975 
rate.  Our  estimates  range  from  -0.13  k to  -0.35  k 
I hc  greenhouse  effect  ol  CFt'li  and  CF  t'l  . on  the 
other  hand,  tends  to  warm  the  surface  by  an  estim- 
ated 0.44  k.  The  net  effect  on  surface  temperature  is 
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Table  1 1 . Factors  used  in  calculating  the  change  in  surface  temperature  due  to  an  SO2  emission  rate  of  3.44  X 
107  kg/yr  in  the  Northern  Hemisphere. 


Source 

Altitude  of 
injection, 
km 

Residence 

time, 

years 

Conversion 

efficiency 

Fraction  in 
hemisphere 

At 

ATs,a 

K 

Pollack  el  al., 

1976 

17 

0.978 

0.869 

0.7 

3.9  X I0-4 

-0.002  to 
-0.004 

At  = 0.031  mb 

20 

1.715 

0.945 

0.7 

7.6  x I04 

-0.005  to 
-0.008 

CTAP  Report  of 
Findings,  1976 

17 

1.70 

0.83 

1.0 

1.1  X I0-3 

-0.07  to 
-0.0  If 

At  = 0.038  m 

20 

3.00 

0.93 

1.0 

2.3  X 10-3 

-0.014  to 
-0.023 

*ATs  = (-6.3to-10K)  at. 

hm  is  (he  mass  density  of  sulfate  aerosols  in  g/m3  averaged  over  a 10-kni  thick  shell. 


uncertain,  but  it  may  be  in  excess  of  ±0. 1 K,  which 
would  be  considered  significant.  The  potential 
climatic  effect  of  a CFM  perturbation  needs  further 
study  to  reduce  the  large  uncertainty  associated  with 
these  results. 


3.6  Effect  of  Ozone  Reductions  on 
Erythema  Dose 

Because  reductions  in  total  ozone  would  permit 
greater  amounts  of  uv  radiation  to  reach  the  surface 
of  the  earth  (Cutchis.  1974;  Halpern  el  al..  1974).  a 
number  of  studies  have  been  performed  with  the  goal 
of  assessing  biological  sensitivity  to  ozone-induced 
changes  in  uv  radiation  (National  Research  Council, 
1973).  One  approach  has  been  to  correlate  skin 
cancer  incidence  data  directly  with  ozone  layer 
thickness.  The  possible  influence  of  such  factors  as 
duration  of  sunlight,  clothing  and  exposure  habits, 
and  optical  path  length  have  been  considered 
(McDonald.  1971;  van  der  Leun  and  Daniels.  1975). 
Another  approach  has  been  to  explicitly  consider  the 
dose  of  uv  radiation  received  as  a function  of  ozone 
amount  and  other  climatic  variables  (Green  and  Mo. 
1975).  The  radiation  dose  is  then  related  to  cancer 
incidence  after  weighting  by  a wavelength-dependent 
function  accounting  for  variation  in  radiation 
efficacy.  This  second  approach,  though  less  direct,  is 
appealing  because  the  mechanism  of  cancer  produc- 
tion is  more  fully  represented,  and  it  allows  for 
experimentation  with  combinations  of  independent 
variables  outside  the  rather  narrow  range  of  reliable 
observation. 


As  the  first  step  in  this  approach.  Green  and  co- 
workers  developed  a semi-empirical  model  for 
calculating  uv  radiation  at  the  surface  of  the  earth  in 
the  spectral  region  280-340  nm  (Green  el  al..  1974a. b; 
Mo  and  Green.  1974).  They  have  calculated  ery- 
thema (sunburn)  dose  as  a function  of  total  ozone, 
solar  zenith  angle,  latitude,  season,  and  cloud 
amount.  Their  model  was  also  used  in  the  Climatic 
Impact  Assessment  Program’s  analysis  of  ozone 
depletion  (Green  el  al.,  1975). 

One  feature  of  this  model,  and  uv  dosimetry  in 
general  (Berger  el  al..  1975).  is  that  the  receiver  is 
assumed  to  be  horizontal.  The  geometry  of  humans, 
however,  is  such  that  the  majority  of  exposed  skin 
would  normally  be  in  a nonhorizontal  position.  In 
fact,  the  horizontal  projection  of  an  upright  person 
amounts  to  only  a few  percent  of  total  surface  area 
(Fanger,  1970). 

The  report  of  the  National  Research  Council 
(1976a)  summarizes  the  specific  sites  of  origin  of 
primary  melanoma  skin  cancer.  The  data  show  the 
predominant  localization  of  sites  oforigintoareasof 
the  skin  that  are  constantly  exposed,  such  as  the 
trunk  and  legs  ol  males.  These  are.  in  general  non- 
horizontal surfaces,  which  suggests  that  receiver 
orientation  should  be  considered  in  the  study  of  uv 
dose. 

In  last  year's  annual  report  (Luther  ei  al..  1977). 
we  reported  the  results  of  a studs  of  the  effect  of 
receiver  orientation  on  erythema  dose.  The  work 
reported  here  expands  upon  that  earlier  effort 
Although  the  same  numerical  model  is  used  for  the 
current  study,  we  will  repeat  the  description  of  the 
model. 
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Model  Description 

I he  procedure  employed  here  involves  using  the 
Cireen  model  (Green  el  al . 1974a)  to  obtain  total 
local  u\  radiation  (Q  + qD  as  a (unction  ol  wave- 
length (Al  and  solar  zenith  angle  (0).  \ ahd  lot  clear 
skv  conditions,  it  is  essentially  a Beer-Bouguer 
tormulation  lor  both  direct  beam  (Q„)  and  ditluse 
(q. ) radiation  It  is  assumed  that  the  latter  is  isotro- 
pically distributed.  I hese  I luxes  are  determined 
using  the  relationships 

y.(0.A>  = H(  A)  \ exp[-A,(0.A)]  cos  0 ( I) 

and 

q.(0.A)  H(  Al  V e.\p(-D,(0.A)]  . (2) 

where  H(A)  is  the  extraterrestrial  solar  irradiance  at 
the  mean  earth-sun  distance  (taken  from  Howard  el 
al.,  I960).  I he  symbol  V.  which  does  not  appear  in 
the  Green  model,  is  the  ratio  of  mean  earth-sun 
distance  to  that  distance  on  a particular  day.  It  thus 
allows  for  eccentricity  in  the  earth's  orbit,  which 
causes  variations  in  the  extraterrestrial  flux  ot 
±3.5rf. 

I he  terms  A and  I),  are  optical  thickness  tunc- 
tions  accounting  for  the  presence  ol  ozone,  air.  and 
particulate  matter,  and  they  are  giv  en  by  Green  <7  a I. 
( 1 974a ).  For  the  spectral  region  considered  here,  the 
total  ozone  optical  depth  is  of  special  interest. 

I o calculate  the  direct  beam  radiation  incident 
on  an  inclined  surface  (Qr).  it  is  necessary  to  replace  0 
in  Fq.  ( I ) w ith  the  angle  between  the  position  of  the 
sun  and  a line  normal  to  the  surface  (i). 

Q,.(0.  A) 

= Him  V exp[-A,(0.A)]  cos  i . (3) 

w here 

cos  i = cos  o cos  0 

+ sin  <»  sin  0 cos  (a  - a')  . (4) 

Here  a'  is  the  solar  azimuth  (relative  to  north),  a is  the 
azimuth  of  the  receiver,  and  a is  the  angle  ol  the 
receiver's  inclination  (relative  to  horizontal).  For  a 
horizontal  receiver,  o 0.  Fquations  required  to  find 
0 and  a’  as  a function  of  latitude,  season,  and  time  ol 
day  may  be  found  in  Dave  el  al.  ( 1975).  Sellers)  1965) 
or  the  Smithsonian  Meteorological  Fables  (Fist. 
195b)  Fhe  last  reference  also  contains  tabled  values 
of  V'. 

Diffuse  radiation  incident  on  an  inclined  surface 
has  two  components:  that  arriving  from  the  sky  as 


scattered  radiation  (q.)  and  that  reflected  front  the 
earth's  surface  (qe).  Under  (he  assumption  of  iso- 
tropic reflection  and  the  earlier  assumption  ol  iso- 
tropic atmospheric  scattering,  these  quantities  are 
given  by  (Dave  <7  al..  1975): 

qs(0.  A)  = 4 qh(0.A)  ( 1 + cos  a)  (5) 

and 

g,.(0.X)  = ^-  |Qh(0,X)  + qh(0,A)|  (1  -cosa).  (6) 

I he  ground  reflectivity  r is  small  for  surfaces  other 
than  snow  and  ice  at  uv  wavelengths  (Kondratyev. 
1973).  A value  of  0.1  was  chosen  here  for  r. 

I he  sum  Q,.  + q,  + q,  represents  the  total  radia- 
tion incident  per  unit  area  on  a plane  inclined  at  an 
angle  a and  with  azimuth  a.  Because  radiation  is  not 
equally  efficient  at  all  wavelengths  in  producing  a 
biologic  response,  it  is  not  possible  to  integrate  inci- 
dent radiation  directly  over  wavelength  and  obtain  a 
meaningful  measure  of  skin  tissue  insult.  It  is  neces- 
sary to  know  the  relative  response  of  skin  to  uv 
radiation  as  a function  of  wavelength.  In  the  case  of 
sunburn  the  so-called  action  spectrum  (or  ery  thema 
efficiency)  has  been  studied  using  uv  lamps  (e  g.. 
C'oblentz  and  Stair.  1934;  Mangus.  1964;  Crippsand 
Ramsey  . 1970)  and  can  be  represented  by  (Green  et 
al..  1974b): 


Fhe  action  spectrum  for  skin  cancer  in  humans  is  not 
known;  it  is  usually  assumed  to  be  the  same  as  for 
erythema  or  similar  to  the  DNA  action  spectrum 
(Setlovv.  1974). 


Procedure  and  Results 

Equations  (4)  through  (7)  may  be  used  to  calcu- 
late the  instantaneous  erythema  dose  for  any  time  of 
day.  date,  latitude,  ozone  amount  and  receiver 
position.  In  analyzing  the  effects  ol  receiver  orienta- 
tion. we  have  numerically  integrated  over  wave- 
length and  time  to  produce  daily  ery  thema  doses  for 
Northern  Hemisphere  latitudes  assuming  the  re- 
ceiver is  stationary . Integrations  were  performed  for 
the  15th  day  of  each  month  using  a time  step  of 
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approximately  20  minutes.  The  total  ozone  was 
specified  as  a function  of  latitude  based  on  data  from 
the  Nimbus  III  satellite  reported  by  Lovill  (1972). 
The  values  of  total  ozone  used  in  the  calculation  are 
shown  in  fable  12.  Ozone  reductions  of  10  and  20% 
were  also  considered  in  the  calculations. 

The  azimuth  of  the  receiver  was  specified  in  two 
ways.  First,  to  obtain  a measure  of  the  uppv  bound 
on  erythema  dose  at  middle  latitudes  in  the  Northern 
Hemisphere,  the  receiver  was  held  fixed  in  a south- 
facing  position.  Second,  to  estimate  the  average  dose 
for  a population  where  there  is  random  orientation 
(i.e.,  no  preferred  orientation),  calculations  were 
made  while  averaging  over  azimuth  angle  (a  rotated 
receiver).  Doses  were  computed  at  20'  intervals  m 
the  azimuth  angle,  and  the  average  of  these  was  used 
for  integration.  The  inclination  angle  n was  held 
fixed  at  values  of  0.  45.  and  90  degrees.  Introducing 
time  dependence  for  this  term  would  require  model- 
ing specific  body  sites  and  activities;  we  have  chosen 
to  present  more  general  calculations. 

Figure  54  shows  the  daily  total  erythema  dose 
for  a south-facing  surface  as  a function  of  latitude  in 
the  Northern  Hemisphere  for  different  ozone 
amounts  and  receiver  inclinations  for  the  1 5th  day  of 
June.  September,  and  December.  I he  relationship 
between  dose  magnitude  and  inclination  is  such  that 
inclined  surfaces  generally  receive  less  daily  inte- 
grated dose  than  a horizontal  surface  at  the  same 
latitude.  Exceptions  to  this  occur  in  some  low  sun 
situations  (large  0)  where  the  dose  on  the  inclined, 
south-facing  surface  slightly  exceeds  that  of  the 
horizontal  surface  (e.g.,  a = 45°,  50-60° N in 
September  and  a = 45°.  20-30°  N in  December).  This 
result  is  in  contrast  to  the  study  of  Dave  el  at.  ( 1975) 
using  the  full  solar  spectrum  where  the  total 
radiation  on  south-facing  surfaces  was  shown  to 
grow  larger  than  that  received  by  a horizontal  surface 
at  high  latitudes.  At  the  wavelengths  contributing  to 
erythema,  the  radiation  is  predominantly  diffuse. 
Consequently,  the  angle  of  incidence  does  not  exert 
as  strong  an  influence  as  it  would  when  the  radiation 
is  predominantly  in  the  direct  beam,  as  was  the  case 
in  the  study  by  Dave  el  al.  (1975). 


Figure  34  also  indicates  that  the  daily  erythema 
does  not  change  uniformly  with  inclination  angle. 
For  example,  changing  a from  0 to  45°  results  in  less 
change  in  the  daily  erythema  dose  than  results  from 
changing  a from  45  to  90°. 

Variations  in  a affect  the  latitudinal  gradient  of 
erythema  dose.  Increasing  a decreases  the  latitudinal 
gradient  in  erythema  dose  poleward  of  the  peak  dose. 
There  is  a small  reduction  in  the  gradient  for  o = 45° 
as  compared  to  a = 0°,  but  the  gradient  is  reduced  by 
approximately  half  for  o = 90°. 

For  a horizontal  surface,  the  latitude  of  the  peak 
daily  erythema  dose  corresponds  closely  to  the  solar 
declination  angle.  As  n increases  there  is  a poleward 
shift  in  the  latitude  of  the  peak  dose.  In  June,  for 
example,  as  « is  changed  from  0 to  90°,  the  latitude  of 
peak  dose  moves  from  20°  N to  approximately  40°  N . 
A shift  of  similar  magnitude  occurs  in  September. 
I his  shift  is  the  result  of  several  factors. 

I he  inclination  angle  affects  the  amount  of 
direct  solar  flux  received  by  the  inclined  surface.  At 
latitudes  toward  the  equator  from  the  latitude  of 
solar  declination,  increasing  n results  in  a decrease  in 
the  direct  solar  flux.  At  higher  latitudes  there  may  be 
either  an  increase  or  a decrease  depending  on  the 
value  of  or. 

As  already  mentioned,  the  contribution  to 
erythema  dose  due  to  diffuse  radiation  has  two  com- 
ponents: that  arriving  from  the  sky  (Eq.  5)  and  that 
reflected  from  the  earth’s  surfaec(Eq.  6).  An  inclined 
surface  receives  less  diffuse  radiation  from  the  sky 
with  increasing  angle  of  inclination  because  less  sky 
is  visible.  For  example,  a vertical  surface  receives  half 
the  flux  of  that  received  by  a horizontal  surface, 
assuming  isotropic  diffuse  radiation.  On  the  other 
hand,  as  a increases  there  is  an  increase  in  the  diffuse 
flux  component  from  the  earth's  surface.  This  llux 
component  is  generally  small,  however,  because  of 
the  low  surface  albedo  assumed. 

Ehe  net  result  of  the  various  factors  just 
discussed  is  that  there  is  a change  in  the  partitioning 
of  the  direct  and  diffuse  (lux  contributions  to 
erythema  dose  as  « is  varied.  This  is  illustrated  lor  a 
south-facing  surface  in  Fig.  35  which  shows  the 


Table  12.  Values  of  total  ozone  inatm-cm. 


latitude  (degrees  north) 


Months 

EQ 

10 

20 

30 

40 

50 

60 

70 

80 

June 

0.256 

0.259 

0.274 

0.295 

0.326 

0.362 

0.391 

0.400 

0.410 

September 

0.249 

0.246 

0.262 

0.295 

0.316 

0.344 

0.364 

0.372 

0.378 

December 

0.256 

0.256 

0.273 

0.303 

0.334 

0.366 

0.393 

0.399 

0.402 
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I l(i  .35.  I he  fraction  of  the  total  daiU  or  > thema  dove  for  a south 
facing  surface  due  to  diffuse  radiation  (sk>  plus  ground  reflection!: 
(a)  June  15.  (h)  September  15.  amt  (c!  December  15. 

traction  ol  the  total  erythema  dose  due  to  the  diffuse 
llu\  (sk\  plus  ground  reflection).  In  low  sun  situa 
lions  (eg.  at  high  latitude),  increasing  n causes  a 
decrease  in  the  diffuse  fraction  ol  the  total  erythema 
dose,  whereas  the  opposite  is  true  in  the  case  of  a high 
sun  (small  l>) 

I'igure  .lb  shows  the  daily  erythema  dose  tot  a 
rotated,  inclined  surface  as  a function  ot  latitude  lot 
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the  >.inu-  values  ol  total  o/one  as  used  in  fig.  14  I he 
dose  received  h\  the  i Dialed  surface  never  exceeds 
that  for  a horizontal  surface,  whereas  u is  possible  ai 
limes  for  a south  facing,  inclined  surface  to  receive  a 
shghtlv  greater  dose  than  the  horizontal  surface 
I here  is  not  a poleward  shtlt  in  the  latitude  of 
the  peak  erythema  dose  with  increasing  n as  was  the 
case  lor  the  south-facing  surface.  Under  the  assump- 
tion of  isotropic  diffuse  radiation,  the  contribution 
to  erythema  dose  due  to  the  total  diffuse  flux  is  the 
same  for  both  the  south-facing  and  rotated  surfaces 
for  a given  n and  latitude  Consequently,  differences 
between  figs  .14  and  lb  are  due  to  differences  in  (he 
direct  flux  component.  I he  direct  flux  component  is 
generallv  largei  for  the  rotated  surface  than  for  the 
south-facing  surface  at  latitudes  near  the  equator, 
and  the  converse  is  true  at  high  latitudes.  In  middle 
latitudes  (.1(1- 5(1"  N)  the  daily  erythema  dose  is 
smaller  lor  the  rotated  surface  than  for  the  south- 
facing  surface. 

I igure  57  shows  the  daily  erythema  dose 
averaged  ovei  a year  for  a south-facing  surface 
assuming  cloudless  conditions.  In  middle  latitudes 
the  daily  average  erythema  dose  on  a surface  w ith  o 
45"  ranges  from  40  to  Wf  of  the  dose  on  a horizontal 
surface  1 he  dose  is  significantly  less  for  inclination 
angles  greater  than  45". 

figure  IS  shows  the  daily  erythema  dose 
averaged  ov  er  a year  lor  a rotated  surface.  In  this  case 
the  daily  average  erythema  dose  in  middle  latitudes 
on  a surface  vv  it  ho  45"  is  approximately  S.Vr  oldie 
dose  on  a horizontal  surface  At  higher  latitudes  the 
dose  on  the  rotated  surface  is  significantly  less  than 
that  on  the  south-facing  surface. 


I Hi.  V7  Dalis  rnthrma  ilosr  airragrU  orrr  a war  for  a south 
faring  wirfarr. 
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A semilog  .scale  is  used  in  figs.  .17  and  7S  to 
demonstrate  the  effect  of  inclination  angle  on  the 
doubling  distance  for  erythema  dose.  In  middle  lati- 
tudes the  computed  erythema  dose  on  a horizontal 
surface  doubles  over  1 1>  ol  latitude  moving  toward 
the  equator.  I he  angle  ol  inclination  has  vety  little 
elfect  on  the  doubling  distance  as  evidenced  In  the 
similar  slopes  ol  the  curves  in  I tgs  77  and  7S 

I lit*  computed  doubling  distance  is  consistent 
w it h measurements  ol  annual-average  erythema  dose 
(Urbach  and  Davies,  1075).  hut  it  differs  trom  the 
doubling  distance  for  the  incidence  ol  'kin  cancer. 
Data  on  the  incidences  ol  skin  cancel  in  middle  l.iti 
(tides  show  a doubling  over  S to  17  of  latitude 
(National  Research  Council,  l‘)7bb),  indicating  that 
there  are  other  factors  contributing  to  the  incidence 
rate.  Consequently,  predicting  the  ettect  ol  a reduc- 
tion in  total  ozone  on  the  skin  cancel  incidence  rate  is 
more  complex  than  pist  assessing  the  change  in 
erythema  dose  and  scaling  proportionately. 

I he  effect  ol  a reduction  in  total  ozone  ol  1 0‘ ( 
ott  the  annual-average  daily  erythema  dose  lot  a 
rotated  surface  is  shown  in  l ig  W I lie  results  were 
almost  identical  lor  the  south  lacing  surlacc  I he 
amplilication  lactoi  on  ei ytliema  dose  (Adose  AIM 
varies  from  1.7  to  I 4 at  low  latitudes,  from  I Mo  7 0 
at  middle  latitudes,  anil  is  approximated  < at  high 
latitudes  It  is  a property  ol  atmospheric  transmis 
sion  that  a I O’  ( reduction  in  the  ozone  column  causes 
the  greatest  percent  change  in  transmission  lot  the 
largest  optical  depth  Because  the  ozone  column  and 
solar  zenith  angle  both  increase  with  latitude,  the 
amplilication  lactoi  is  laigest  at  high  latitudes  \\  tide 
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l lti.  W.  I hi*  |>rrcrnt  change  in  ilailv  ervtheiiia  dusr  Inr  a .maleil 
surface  averaged  uvrr  a year  ilur  In  a 10'  i leducliun  in  nfiinr. 


the  percent  increase  in  daily  erythema  dose  ai  high 
latitudes  is  much  greater,  the  total  amount  ot  radia- 
tion is  very  small  compared  to  lower  latitudes  Since 
the  maionly  ol  the  world's  population  lives  south  ol 
ss  \ . amplilication  lactoi s in  the  lange  I Mo  2.0  are 
most  realistic 

Reducing  the  ozone  column  incicases  both  the 
direct  and  diltuse  flux  components,  but  by  different 
luclois  Since  the  relative  contribution  ol  the  direct 
and  diltuse  flux  components  to  erythema  dose 
depends  upon  a.  the  amplification  factor  also 
depends  upon  o I his  dependence  is  lessened  at  high 
latitudes  where  the  flux  is  almost  entirely  diffuse 

As  mentioned  earlier,  the  action  spectrum  for 
skin  cancel  is  an  unknown  tunclion.  However,  van 
dei  I cun  and  Daniels  (1075)  argue  that  the  peak  in 
the  action  spectrum  should  occur  neat  770  mu.  I his 
is  slutted  consideiablv  in  wavelength  from  the  loca- 
tion ol  the  peak  in  the  eivthema  action  spectrum 
described  bv  I q (7),  which  occurs  at  297  nm.  In 
ordei  to  test  the  sensitivity  ol  dose  to  action  spec- 
trum position,  the  preceding  calculations  were 
repeated  substituting  770  nm  in  l q.  (7)  Aside  from 
the  expected  decieaxe  in  dose  magnitude,  very  little 
changed  (torn  the  icsiills  obtained.  An  important 
exception  concerns  the  percent  change  in  the  annual 
aveiage  daily  eivthema  dose  pioduced  by  a 10', 
reduction  in  ozone,  which  is  shown  in  f ig  40  Dose 
amplilication  lactois  ate  sigmltcantlv  liighei  at 
middle  and  low  latitudes  computed  to  the  results  m 
I ig  W I he  amplilication  factor  langes  tiom  I S to 
1 0 at  low  latitudes  and  liont  I 0 to  7 t at  middle 
latitudes  I he  reason  lot  the  target  amplification 
lactois  is  related  to  the  ozone  absorption  cross 


I K.  40.  The  percent  change  in  dail)  erythema  dose  for  a rotated 
surface  aseraged  oser  a sear  seith  the  peak  in  the  action  spectrum 
shifted  to  270  nm. 


section,  which  in  the  region  near  300  nm  is  a rapidly 
decreasing  function  of  wavelength.  Shifting  the 
action  spectrum  to  shorter  wavelengths,  where  the  uv 
fluxes  are  more  sensitive  to  variations  in  o/one 
amount,  leads  to  larger  amplification  factors. 

Figure  41  shows  the  distance  moved  south, 
which  is  equivalent  to  a 1 0ri  reduction  in  o/one 
based  on  the  annual  average  data  for  a rotated  sur- 
face shown  in  Fig.  38.  Fhe  shaded  area  indicates  the 
range  of  values  as  a is  varied  fromO  to  90°,  the  largest 
values  corresponding  to  a = 90° . In  middle  latitudes 
the  increase  in  daily  average  erythema  dose  due  to  a 
I O'-;  o/one  reduction  is  roughly  equivalent  to  moving 
south  a distance  of  350  to  450  km  with  no  o/one 
perturbation.  The  equivalent  distance  is  much 
greater  at  latitudes  toward  the  equator  from  30°N. 

I here  is  no  equivalent  distance  for  latitudes  less  than 
I5°N.  because  the  erythema  dose  at  these  latitudes 
with  a I0'r  o/one  reduction  isgreater  than  the  unper- 
turbed erythema  dose  at  the  equator. 


Kl<;.  41.  the  distance  muted  south  equitalenf  to  a IlHi  o/one 
reduction  computed  for  a rotated  surface.  The  shaded  area  indi- 
cates the  range  of  talues  as  the  inclination  anglr  is  taried  from 
0 to  90  . 

Summary 

According  to  the  simplified  model  used  here, 
erythema  doses  have  a significant  dependence  upon 
receiver  orientation.  Although  it  is  possible  at  times 
for  an  inclined  surface  to  receive  a larger  erythema 
dose  than  a horizontal  surface,  in  general  the  highest 
doses  are  received  by  horizontal  surfaces.  In  middle 
and  high  latitudes,  inclined  surfaces  (south-facing  or 
rotated)  with  inclination  angles  up  to  45°  receive  at 
least  80c;  of  the  daily  erythema  dose  received  by  a 
horizontal  surface.  Much  larger  reductions  in  daily 
erythema  dose  occur  for  a > 45°.  which  can  result  in 
a significant  reduction  in  the  latitudinal  gradient  of 
erythema  dose.  In  the  case  of  a south-facing  surface, 
the  inclination  angle  affects  the  latitude  of  the  peak 
erythema  dose.  Varying  the  inclination  angle, 
however,  does  not  have  a significant  effect  on  the 
doubling  distance  for  average  daily  erythema  dose. 
Finally,  if  the  skin  cancer  action  spectrum  peaks  near 
270  nm  as  argued  by  van  der  Leun  and  Daniels 
( 1975).  skin  cancer  rates  may  be  more  sensitive  than 
erythema  to  variations  in  total  o/one. 


4. 


WORK  IN  PROGRESS 


We  are  developing  a new  solution  technique  for 
the  one-dimensional  model.  The  new  code  is  being 
designed  so  that  certain  species  or  species  families 
can  be  held  constant  while  other  short-lived  species 
arc  calculated.  We  may  then  use  the  program  to 
study  how  latitudinal  variations  in  species  concentra- 
tions and  temperature  feedback  control  OH  concen- 
trations in  the  troposphere.  We  expect  to  use  this 
program  to  test  the  sensitivity  of  global  budget 
analyses  to  latitudinal  variations. 

The  new  solution  method  is  also  being  designed 
to  allow  either  steady-state  or  time-dependent 
calculations.  It  may  be  possible  to  extend  the  method 
to  the  two-dimensional  model.  The  basic  solution 
scheme  is  as  follows  (this  is  subject  to  change,  since 
the  method  has  not  been  tested): 

1.  Predict  values  of  individual  species  con- 
centrations at  next  time  step,  t + At. 

2.  Predict  values  of  species  family  concentra- 
tions (e.g.,  NO,.  CIO,,  etc.). 

.V  Use  the  concentrations  from  ( I ) to  estimate 
the  time  derivative  and  llux  terms  of  the 
species  conservation  equation.  Scale  each 
specie’s  time  derivative  by  the  predicted 
time  derivative  for  its  family  from  (2).  This 


will  damp  fluctuations  in  short-lived  species 
concentrations  and  allow  a full  solution  of 
the  stiff  system  of  equations  we  wish  to 
solve.  As  a fully  consistent  solution  is  ap- 
proached. the  scaling  factor  becomes  equal 
to  unity  and  the  true  solution  will  be 
approached. 

4.  Solve  the  system  of  linear  algebraic  equa- 
tions that  result  from  using  the  fixed  terms 
estimated  in  (3)  and  the  individual  produc- 
tion and  loss  terms  for  each  species.  The 
result  is  a new  set  of  species  concentrations 
at  time  t + At. 

5.  Check  to  see  if  the  concentrations  from  (4) 
are  the  same  as  those  from  ( 1 ).  If  not.  use  the 
new  values  from  (4)  and  repeat  the  proce- 
dure starting  with  step  (2). 

A special  code  is  being  written  that  will  enable 
the  model's  vertical  grid  structure  to  be  easily 
changed.  The  first  utilization  of  this  code  will  be  to 
develop  a more  detailed  model  of  atmospheric  chem- 
ical and  dynamical  processes  in  the  troposphere.  This 
model  will  be  used  to  derive  a global  distribution  of 
OH  to  be  compared  with  recent  measurements  by 
Poug  Davis. 
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APPENDIX  A.  CHANGES  TO  THE  ONE-DIMENSIONAL 
TRANSPORT-KINETICS  MODEL  DURING  THE  PAST  YEAR 


The  chemical  reactions  incorporated  in  our  cur- 
rent model  are  given  in  Table  A- 1 along  with  the  rate 
expressions.  Nearly  all  of  the  rate  expressions  come 
from  Hampson  and  Garvin  (1978).  Most  of  the 
exceptions  reflect  recently  reported  measurements, 
although  a few  minor  reactions  not  treated  by 
Hampson  and  Ciarvin  (1978)  come  from  other 
compendia. 


We  have  adopted  the  measurements  of  Howard 
el  al.  (1 978)  for  the  reactions 

O,  + HO;  - OH  + 20; 

and 

NO  + HO;  - NO;  + OH 


Table  A-l.  Chemical  reactions  and  rate  coefficients  currently  used  in 

the  one-dimensional  transport-kinetics  model. 

Reaction 

A 

Rate  coefficient  (k  = Ae°  1 ) 

B 

Ref. 

O + ♦ M = O3  ♦ M 

1.1  x to"-'4 

510 

1 

O + O3  = 202 

1.9  X I0"11 

-2300 

l 

O3  ♦ NO  = NOj  ♦ Oj 

2.1  X 10'12 

-1450 

1 

O + NO 2 = NO  * O2 

9.1  X 10'12 

0 

1 

N2O  ♦ Ot*DI  = N2  ♦ 02 

5.5  X 10"" 

0 

1 

N2O  ♦ 0(1D)  = 2NO 

5.5  X 10"1 1 

0 

t 

N ♦ O2  = NO  ♦ O 

5.5  X I0"12 

-3220 

1 

N + NO  = N2  + 0 

8.2  X 10"" 

-lit) 

1 

O('o)  ♦ H20  = 20H 

2.3  X tO"10 

0 

1 

Of*D)  ♦ CH4  * OH  ♦ 2HOj  ♦ CO* 

1.3  x It)-10 

0 

1 

O3  + OH  = H02  ♦ O2 

1.5  X I0"12 

-1000 

1 

O ♦ OH  = O2  ♦ H 

4.2  X It)"11 

0 

1 

O3  ♦ HO2  = OH  ♦ 202 

1.4  X to"14 

-590 

2 

0 ♦ HO2  = OH  + O2 

3.5  X tO'11 

0 

1 

H *02  ♦ M * Ht>2  ♦ M 

2.08  x I0~i2 

290 

1 

O3  ♦ H = OH  ♦ O2 

1.42  X 10"10 

-478 

H02  ♦ HOi  = H2t)2  * 02 

Note  1 

4 

IK>2  ♦ OH  = H20  ♦ O2 

3.0  X I0"11 

0 

1 

OH  ♦ NO2  ♦ M = HNO3  + M 

Note  2 

1 

OH  ♦ HNO3  = H2O  + NO3 

8.9  x It)"14 

0 

5 

H2O2  ♦OH  = H2  ♦ HO2 

1.0  x to-11 

-750 

1 

Nj  ♦ 0(*D)  ♦ M * N20  + M 

3.5  X tO-17 

0 

1 

N ♦ NO2  = N2O  ♦ O 

2.0  X I0"11 

-800 

1 

NO  + O ♦ M = NO2  ♦ M 

1.6  X tO'42 

584 

1 

NO  ♦ HO 2 = N02  ♦ OH 

3.3  X 10"12 

254 

2 

H2  ♦O('O)  = 0H  ♦ H 

9.9  X I0"11 

0 

1 

OH  ♦ OH  = H2O  ♦ O 

1.0  X to"11 

-550 

t 

N ♦ O3  ■ NO  ♦ 02 

2.0  X 10"11 

-1070 

1 

NO 2 +03 « NO3 ♦ 02 

1.2  X I0"1-’ 

-2450 

1 

OH  ♦ CH4  = H2t)  ♦ 21102  + 00* 

2.36  X I0“12 

-1710 

1 

OH  ♦ OH  ♦ M = Hj02  + M 

1.2  X I0~i2 

900 

1 
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Table  A-l.  (Continued) 

Kite  coefficient  (k  = Ae°  * ) 

K nction 

A 

B 

Ref. 

HjOj  ♦ O = OH  ♦ HO j 

2. 75  X 10"' 2 

-2125 

1 

O ♦ l II4  = OH  ♦ 2 HO  j ♦ 1X1* 

3.5  X I0"M 

-4550 

1 

Cl)  ♦ OH  = H ♦ CX>2 

Note  3 

6 

ot'm  ♦ M *0  ♦ M 

2.2  x 10'" 

*»»» 

1 

O ♦ O3  = CIO  ♦ 02 

2.7  x 10'" 

-257 

1 

O ♦ OOO  ♦ 200 

5.*»  X It)'" 

0 

1 

O ♦ CH4  = HO  ♦ 2H02  ♦ CX)1 

7.3  x I0'12 

-1260 

1 

O ♦ NO2  ♦ M = CINO2  ♦ M 

ti.'i  X I0'M 

2115 

7 

O + CINO2  = 20  + NO2 

3.0  x It)"12 

0 

1 

00  ♦ O = Cl  ♦ O2 

7.7  x It)"1 1 

-130 

1 

NO  ♦ OO  S NO2  + Cl 

1.0  x 10'" 

200 

1 

00  ♦ O3  = a ♦ 02  + 02 

1.0  x I0'12 

-4000 

1 

OO  ♦ O3  = OOO  ♦ 02 

1.0  X I0"12 

—4000 

1 

OO  ♦ N<>2  ♦ M = OONO2  ♦ M 

Nolo  4 

1 

OO  ♦ CK>  = O + OCR) 

2.1  x I0~12 

-2  200 

1 

OO  + OO  = 20  ♦ l>2 

1.5  X I0"12 

-1238 

1 

HO  + Ot'l))  » O + OH 

1.4  x I0'10 

0 

1 

OH  ♦ HO  = 112  ♦ 0 

3.0  X I0'12 

-425 

1 

0 ♦ HO  = OH  ♦ 0 

1.14  x It)"1 1 

-3370 

1 

O ♦ OOO  = OO  ♦ l>2 

2.0  x 10'" 

- 1 1 00 

1 

NO  ♦ OCIO  = Nt>2  ♦ OO 

3.4  x I0"n 

0 

8 

N * OCIO  a no  + CK) 

K.O  x It)"1  i 

0 

8 

H ♦ OOO  = OH  ♦ CK) 

5.7  X 10'" 

0 

1 

O ♦ OH  = HO  ♦ O 

1.0  x 10'" 

-2470 

1 

O ♦ IU>2  = HCI  ♦ 02 

4.1  x 10'" 

0 

4 

O ♦ HNO3  = HCI  ♦ NO3 

1.0  x It)'1 1 

-2170 

1 

CFCI3  ♦ O('d)  = 2Clb 

2.3  x I0~10 

0 

1 

CF2CI2  + 0(lD)  = 2Clb 

2.0  x 10-'° 

0 

1 

tl  + 112  = HCI  ♦ H 

4.7  x 10'" 

-2340 

It) 

a ♦ H2O2  * HCI  ♦ HO2 

1.6  x I0"12 

-384 

1 

OONO2  ♦ O = CK)  ♦ NO3 

3.0  x It)'12 

-808 

1 

CH3CI  ♦ OH  = Cl  ♦ II2O  ♦ tU)2b 

2.2  x I0'12 

-1142 

1 

NO  ♦ NO3  = 2NO2 

8.7  X It)'12 

0 

II 

NO2  ♦ 0 + M = NO3  ♦ M 

1.0  x 10'-" 

0 

1 

NO2  ♦ NO 3 = NO  ♦ O2  + NO2 

2.0  x 10',J 

- 2000 

1 

NO 2 ♦ NO3  + M = N2O5  ♦ M 

Note  5 

6 

N2O5  ♦ M = Nt)2  + NO3  ♦ M 

Nolo  t 

6 

N2O5  ♦ O = 2NO2  ♦ Oj 

1.0  x I0'14 

0 

12 

n2°5  ♦ H2O  = 2HNO3 

1.0  x It)"20 

0 

1 

Oi'd)  ♦ O3  = 202 

1.2  x I0~10 

0 

1 

IK)2  ♦ IU>2  ♦ H2O  = H2O2  + O2  ♦ H2O 

Note  7 

4 

1 


r. 

I 

l 


4(> 


Table  A-l.  (Continued  I 


Photolysis  tractions 

02  = 0 + 0 
Oj  = O ♦ o2 
O3  = 0('d)  + 02 
NO2  = NO  ♦ O 
N2O  = N2  + 0('d> 

NO  = N + O 
HNO3  = OH  ♦ NO2 
H2O2 = 20H 
»R>2  = OH  + O 
CIONOj  = OO  + NO2 
HCI  = H + ci 

cio  = a + 0 
cio  = a + o(,D) 

ON02  = o ♦ NO2 
000  = 00  + 0 

CFjCli  = 2Clb 
CFO  3 = 3Clb 
CCI4  = 4Clb 

CH3CI  = 2H02  + CO  + Ob 
N2O5  =2N02  +0 
NO3  = NO  + O2 
NO3  = NO2  + O 
H2O  = H + OH 


Note  I: 


3.25  v 108  + M X 4 X 10-10 


(I  + 3.5  x 10-16  x Mxe~20b0T)  (8M  + 4 x 10+20) 


Note  2: 
Note  3: 


l s»nji  (1977)  as  given  by  Hudson  (1977). 
1.4  X 10~'J  + Mx  7.3  X I0"J\ 


Note  4: 


3.3  x IO~23  t~3-34 
1.  + 8.7  X 10-9  T~°  6 X (M)1  2 


Note  5: 


2.9  X 10 


12 


7.0  X I0+21  e”2670T  + M 


Note  6: 


6 X 10+**  e~l0970fT 
7.0  X I021  e-267(,rr  + M 


Note  7: 


, , x ,0-34e+3730rr 
1 ♦ M x 3.S  x I0~,be~20b0  1 


*Methyl  radical  production  is  assumed  to  yield  21102  * CO. 

b Products  such  as  F,  H2.  CO.  OOF2.  etc.  are  not  followed  when  prod uced  from  lulocarbons. 
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We  have  adopted  the  pressure,  temperature  and 
water-vapor-dependent  expression  for 

M.H20 

HOs  + H02 ► H202  + 02 

as  given  by  Cox  (1978)  in  his  errata.  We  have  adopted 
values  given  by  Watson  (1978)  for 

H + O,  - OH  + O; 

and  Thrush  (1978)  for 

Cl  + HO;  - HC1  + O;  . 

A privately  communicated  estimate  by  Johnston 
(1977)  is  used  for 

N.O,  + O - 2NO;  + O;  . 

Several  minor  reactions  involving  OCIO  were 
retained  based  on  expressions  cited  by  Watson 
( 1975).  None  of  these  has  any  effect  on  model  predic- 
tions. Simple  pressure-dependent  rate  expressions 
for 


Table  A-2.  Chemical  rate  coefficients  used  in  the  1977 


model  chemistry.  Only  those 
changed  during  the  past  year  are 

reactions  that 
listed . 

were 

Rate  coefficient 
(k  = AeBAF) 

Reaction 

A 

B 

O3  + HO2  -*  OH  + 202 

1.0  X 10'13 

-1 250 

0 + HO2  - OH  + O2 

3.0  X 10-11 

0 

O3  + H-OH+O2 

1 .23  X I0-10 

-562 

HO2  + HO2  - H2O2  + O2 

1.7  x 10-11 

-500 

HO2  + OH-H2O  + O2 

2.0  x I0",J 

0 

NO  + HO2-  NO2  + OH 

4.28  x 10~" 

-500 

N +O3-*  NO  + O2 

4.0  x 10-12 

-650 

CO  + OH  -*  H + CO2 

1.4  X 10“13 

0 

NO  + CIO  - NO2  + Cl 

2.2  x 10"11 

0 

CIO  + NO2  + M -»  CIONO2  + M 

5.1  X 10'33 

1030 

Q + HO2  - HCI  + O2 

3.0  x 10_1 1 

0 

M 

NO2  + NO3-N2O5  , 

M 

N2O5  -*■  NO2  + NO3 

HO  + CO  - C02  + H , 

were  developed  from  data  cited  by  Hampson  and 
Garvin  (1978). 

It  should  be  noted  that  we  do  not  use  a detailed 
mechanism  for  the  oxidation  of  methyl  radicals. 
Instead  we  assume  that  a methyl  radical  is  oxidized  in 
such  a way  as  to  releast  two  HO;  radicals  and  CO. 

Rate  coefficients  were  changed  for  several  reac- 
tions during  the  past  year.  These  reactions  are  listed 
in  T able  A-2  along  with  the  rate  coefficients  that  were 
used  previously.  The  most  significant  changes  are  the 
new  rate  for  HO;  + O3,  the  new  temperature  depen- 
dence for  HO2  + NO,  and  a pressure-  and  water- 
vapor  dependent  rate  for  HO;  + H02. 

Water-vapor  concentration  is  now  computed  at 
altitudes  greater  than  1 2 km  in  the  model.  The  water- 
vapor  distribution  in  the  troposphere  (0-12  km)  is 
specified  based  on  the  global  average  distribution 
given  by  Oort  and  Rasmusson  (1971).  The  derived 
ambient  water-vapor  distribution  increases  in  con- 
centration from  3.5  ppmv  at  the  tropopause  to  5.4 


ppmv  in  the  upper  stratosphere.  Water  vapor  is 
produced  in  the  upper  stratosphere  by  the  reactions 
OH  + H0;-H;0  + 0;and0H+CH4-H;0  + CH,. 

A new  scheme  was  developed  to  calculate  the 
multiple  scattering  factors  used  in  the  calculation  of 
photodissociation  rates.  A computationally  fast 
method  for  calculating  the  multiple  scattering  factor 
was  needed  for  the  diurnal  model,  and  this  method  is 
also  used  in  the  diurnal-averaged  model.  Multiple 
scattering  factors  are  no  longer  read  in  but  are  com- 
puted internally  for  the  given  solar  zenith  angle  and 
species  distributions.  A detailed  description  of  the 
method  will  be  prepared  in  the  near  future. 

An  effort  has  been  made  to  program  the  model 
to  run  as  efficiently  as  possible,  to  use  a minimum 
amount  of  storage,  and  to  exist  in  a modular  form  for 
ease  of  future  modifications.  To  minimize  the  use  of 
storage,  the  model  was  divided  into  sections  with  the 
use  of  overlays  so  that  at  runtime  all  coding  used  to 
initialized  the  problem  is  released.  There  is  now  more 
space  to  add  future  chemistry.  A significant  reduc- 
tion in  running  time  was  accomplished  by  using 
machine  language  versions  of  two  subroutines.  The 
user  now  has  the  option  of  running  either  diurnal  or 
diurnal-averaged  versions  of  the  model  with  or  with- 
out including  the  effect  of  temperature  feedback. 
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